Superconductive proximity effect in interacting disordered conductors 
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We present a general theory of the superconductive proximity effect in disordered normal- 
superconducting (N-S) structures, based on the recently developed |jj Keldysh action approach. 
In the case of the absence of interaction in the normal conductor we reproduce known results for 
the Andreev conductance Ga at arbitrary relation between the interface resistance Rt and the dif- 
fusive resistance Rd- In two-dimensional N-S systems, electron-electron interaction in the Cooper 
channel of normal conductor is shown to strongly affect the value of Ga as well as its dependence on 
temperature, voltage and magnetic field. In particular, an unusual maximum of Ga as a function 
of temperature and/or magnetic field is predicted for some range of parameters Rd and Rt- The 
Keldysh action approach makes it possible to calculate the full statistics of charge transfer in such 
structures. As an application of this method, we calculate the noise power of an N-S contact as a 
function of voltage, temperature, magnetic field and frequency for arbitrary Cooper repulsion in the 
N metal and arbitrary values of the ratio Rd / Rt ■ 

PACS numbers: 74.40. +k, 74.50. +r, 72.10.Bg 



I. INTRODUCTION 

The essence of the superconductive proximity effect 
is an existence of Cooper correlations between electrons 
in the normal metal in contact with a superconductor. 
The microscopic mechanism leading to the proximity ef- 
fect is Andreev reflection M of an electron into a hole 
at a normal-metal-superconducting (N-S) interface. The 
probability of Andreev reflection (as opposed to the nor- 
mal reflection) and thus the strength of the proximity 
effect is determined by the transparency of the N-S in- 
terface, which may be relatively weak due to the presence 
of a tunnel barrier or mismatch between Fermi velocities 
of two materials. Disorder in the normal conductor near 
an N-S contact is shown theoretically (^-j^] to increase 
considerably the effective probability of Andreev reflec- 
tion (see Ref. |?J for a recent review from the experimen- 
tal viewpoint). However, when the normal conducting 
region is made of a dirty metal film, or two-dimensional 
electron gas with low density of electrons, Coulomb in- 
teraction in the normal region (neglected in -|(|) gets 
enhanced [pi and leads to strong quantum fluctuations 
which suppress the Andreev conductance. Several differ- 
ent kinds of quantum effects are known to be relevant in 
low-dimensional conductors at low temperatures. Quan- 
tum corrections to the conductivity of two-dimensional 
dirty conductors grow logarithmically as temperature T 
decreases and become large at ln(l/Tr) ~ g (where 
g = (H/e 2 )a is the dimensionless conductance, a is the 
conductance per square, and r is the elastic scattering 
time) . There are two main types of these effects: weak lo- 
calization corrections , and interaction- induced cor- 
rections Other important quantum effects include 
interaction-induced suppression of the tunneling conduc- 
tance (Coulomb zero-bias anomaly PO]), and disorder- 



induced suppression |12|-|lq| of the superconductive tran- 
sition temperature T c . The corresponding corrections are 
of the relative order of g -1 ln 2 (l/Tr), i. e., much stronger 
than the weak localization and interaction corrections. 

In the previous paper |jj we have studied the influ- 
ence of the last two effects on the Andreev conduc- 
tance and the Josephson proximity coupling in S-I-N and 
S-I-N-I-S structures (where T denotes an insulating tun- 
nel barrier). Using the renormalization group method, 
we have calculated both quantities including the region 
of strong fluctuations, ln 2 (l/Tr) > g. Cooper-channel 
repulsion was found to modify results of semiclassical 
calculations by the power-law factors oc T 1 ^^. The 
effects of the Coulomb zero-bias anomaly (in the case of 
unscreened static Coulomb interaction) were found to be 
even stronger, the corresponding correction factor being 
of the order of exp [-(l/47r 2 ff ) ln 2 (A/T)], where A is the 
gap in the superconductive terminal. However, results of 
Ref. [0 are limited to the lowest tunneling approxima- 
tion, i. e., valid under the condition Rt 3> Rd, where 
Rt and Rd are the resistances of the tunnel barrier in 
the normal state, and of the diffusive normal conductor, 
correspondingly. General semiclassical theory of N-I-S 
systems with arbitrary ratio Rt/ Rd usually neglects in- 
teraction effects in the N part of the structure 
The effect of Cooper interaction (with strength A) on 
the Andreev conductance Ga was studied by Stoof and 
Nazarov for an N-S structure with the ideal interface 
and ID geometry of the normal wire, in the first order 
over A. Relative correction to Ga was found to be small, 
of the order of A <C 1. The effects of interaction upon 
noise in N-S structures had never been studied before, 
to the best of our knowledge (for the discussion of N-S 
noise in the absence of interaction, see Refs. |l8|,[l1| and 
the review pQ]). 
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FIG. 1. A small superconductive island (S) of size 2d con- 
nected to a reservoir (R) through a tunnel barrier (I) and a 
dirty normal film (N) of size 2L 3> 2d. 

In this paper we generalize our approach Jlj for the case 
of strong proximity effect (arbitrary relation between Rt 
and i?£>), in the presence of Cooper interactions in the 
normal conductor. We calculate the Andreev conduc- 
tance of 2D N-I-S structures (shown in Fig. [[]) at low 
(compared to A) temperature and voltages, as a func- 
tion of t = Rd/Rt and of the "decoherence time" of 
an electron and the Andreev- reflected hole, fi./f2*, where 
£1* = ma,x(T,eV,eDH/c). We demonstrate that the 
Cooper interaction effects in 2D are substantially differ- 
ent from the ID case considered in Ref. J^|. In par- 
ticular, at 0* <C Et\i = hD/L 2 the lowest-order relative 
correction 5Ga/Ga scales as A In L/d and grows with the 
size of the system L. We sum up the main logarithmic 
terms of the order of (A In L/d) n and find that Cooper 
repulsion may lead to nonmonotonous dependence of G a 
on Rd/Rt and/or on the decoherence energy scale f2*. 

Technically, our method is based on the Keldysh func- 
tional approach for disordered superconductors jlj] (pro- 
posed previously for normal conductors in Rcf. fell ). 
The basic object of this theory is the 8x8 matrix 
field Q(t,t';r) which depends on two time variables and 
one space coordinate. Its average value gives the elec- 
tron Green function G(r, r) at coincident points. Fluc- 
tuations of the Q matrix describe slow diffuson and 
Cooperon modes of the electron system. In the case of 
spin-independent interactions, Q reduces to a 4 x 4 ma- 
trix. This method seems to be more convenient than 
the replica functional approach [l4| - [l^1 as it does not 
require tedious analytic continuation procedure and al- 
lows for the direct study of non-equilibrium phenom- 
ena. In the limit Rt S> Rd studied in 0], integra- 
tion over diffuson/ Cooperon modes reduces the prob- 
lem to an effective action S plOK [Qs, Qn] for the prox- 
imity effect, which contains two terms only, Tr QsQn 
and Tt(QsQn) 2 , describing transport of single electrons 
and Cooper pairs, respectively (here Qs and Qn corre- 
sponds to the superconductive and normal boundaries 
of the system). At arbitrary relation between Rt and 
Rd, the proximity action contains an infinite series of 
terms, S piox [Q s , Qn] oc Y^Li 7n ^(QsQn) 71 - A set 



of parameters 7„ can be conveniently parametrized via 
the 27r-periodic function u(x) = ~Y^ =1 n r y n smnx that 
encodes amplitudes of processes of n-electron transfer 
through the system. The evolution of u(x) as a func- 
tion of the ratio Rt/Rd is found with the use of the 
functional renormalization group method applied to the 
action S prox [Q s , Qn}- 

In this paper we will not consider the weak localization 
and interaction corrections assuming that the sheet con- 
ductance is relatively large, g ^> \n(L/d). We will also 
not take into account the effect of the Coulomb zero- 
bias anomaly (ZBA) on the Andreev conductance. Pos- 
sibility of neglecting Coulomb ZBA effects while treat- 
ing the Coulomb-induced repulsion in the Cooper chan- 
nel |^4]-|l6|] is due to the fact that the specific form of 
the Coulomb ZBA depends crucially on the long-range 
behavior of the Coulomb potential, whereas renormal- 
ization of the Cooper-channel interaction depends on the 
short-distance Coulomb amplitude only. If long-range 
Coulomb forces are suppressed (i. e., by placing a nearby 
screening metal gate), the Coulomb ZBA effect may be- 
come weak, and the main effect comes from the short- 
range repulsion in the Cooper channel. This is the situ- 
ation we are going to study in this paper. 

Another limitation of the present discussion is that we 
will consider the case of a two-dimensional (2D) geome- 
try of the current flow between the superconductive and 
normal electrodes of the structure, as shown in Fig. [l]. 
This will make possible to construct a unified functional 
renormalization group treatment that takes into account 
modifications of the proximity effect strength both due to 
multiple Andreev reflections and due to Cooper-channel 
repulsion. 

Finally, we restrict our present discussion to a low- 
frequency domain u> <C D/L 2 . It can be shown that the 
frequency-dependent Andreev conductance Ga(w) does 
not contain significant frequency dependence at scales 
u> > D/L 2 , due to the long-range nature of the Coulomb 
interaction which makes electron liquid effectively incom- 
pressible up to much high frequency scale defined by the 
inverse time of electric field propagation across the struc- 
ture. However, the specific form of the proximity ac- 
tion iSproJQs, Qn] as a linear combination of the traces 
of (QsQn)™ is valid only in the low-frequency domain. 
General discussion of the action becomes much more in- 
volved in the high frequency region, and we will postpone 
this subject for the future studies. 

The rest of the paper is organized as follows. In Sec. || 
we introduce the Keldysh cr-model action and describe 
a general scheme of expansion in terms of diffusive slow 
modes. Sec. [II is devoted to the derivation of the func- 



tional renormalization group (FRG) equation for the case 
of non-interacting electrons in the normal conductor. It 
is shown that the resulting equation for the function u(x) 
which parametrizes the action exactly coincides with the 
Euler equation for generating function of transmission 



2 



eigenvalues Tj discussed in Ref. [E2|. In Sec. [V we show 



how physical quantities such as conductivity and current 
noise can be expressed in terms of the function u(x). The 
main new results of the paper are presented in Sec. 
In Sees. VA and VB we derive an additional term in 



the FRG equation, that accounts the electron-electron 
interaction in the Cooper channel. The solution of the 
full FRG equation and the corresponding results for the 
Andreev conductance are discussed in Sec. VC. The ef- 



fect of interaction upon the current noise is studied in 
Sec. VD. Sec. VI is devoted to the study of the de- 



pendence of the conductance and noise on temperature, 
voltage and magnetic field, which suppress Cooperon am- 
plitude at the e nergy scale fi* = max(T, eV, eDH/c). It 
is shown in Sec. VIC that repulsion in the Cooper chan- 
nel may lead to a nonmonotonous dependence of effective 
interface resistance on fi*. Sec. |VII contains discussion 
and conclusions. Some technical details are presented in 
Appendixes |X| and Appendix ^| contains an analysis 
of the Andreev conductance in a 2D S-I-N structure by 
means of the standard method of the Usadel equation in 
the presence of Cooper interaction. 



II. KELDYSH cr-MODEL ACTION 

Keldysh action for disordered N-S systems derived in 
Ref. p| can be represented as a sum of the bulk and 
boundary contributions. 

The bulk action, SWik, is a functional of three fluctu- 
ating fields: the matter field Q(r,t,t') in the film, the 
electromagnetic potential 4>(v,t) and the order parame- 
ter field A(r, t) used to decouple the quartic interaction 
vertex in the Cooper channel: 



5hiilk = 



17TIS 



4 

-Tr~. 



■Tr 



£>(VQ) 2 
T (V ~ 1 + 2v)a x % 



-A)Q 
•A-«x x A : (2.1) 



+ 4i(iT z d t 

2v . 
T 

where D is the diffusion coefficient in the film, v is the 
density of states at the Fermi level per single projection 
of spin, and A is the dimensionless coupling constant in 
the Cooper channel. 

Being a matrix in the time domain, Q(r, t, t') is also a 
4x4 matrix in the direct product K ® N of the Keldysh 
and Nambu spaces. Pauli matrices in the K and N spaces 
are denoted by Oi and Tj respectively. The field Q satisfies 
a nonlinear constraint Q 2 = 1 and can be parametrized 
as Q = U~ l MJ where A = A t z is the metallic saddle 
point and 



A (e) = 



1 2F(e) 
-1 



(2.2) 



K 



The matrix F(e) introduced in Eq. ( |2.2| ) acts in the 
Nambu space and has the meaning of a generalized dis- 
tribution function. In equilibrium F(e) = tanh(e/2T), 



and its general form is F(e) = /(e) + fi(e)r z , where 
fi(E) is the anomalous distribution function which is 
a measure of the branch imbalance |2^,|2^]. The object 
4> = {4>i,4>2) T is a vector in the Keldysh space, with 
4>i , (f>2 being the classical and quantum components of 
the field. They are given by the symmetric and antisym- 
metric linear combination of the <j> fields on the forward 
(f) and backward (b) branches of the Keldysh contour: 
0i,2 = (</>f ± </>t>)/2. In the Keldysh technique, exter- 
nal (nonfluatuating) fields (e. g., applied voltage) have 
only classical component, while physical results can be 
obtained by taking derivatives with respect to the quan- 
tum component, cf. Sec. IV. Fluctuating fields have 
both components, 
matrix (f> — (fiiao 



is a shorthand notation for the 
t>2&x in the Keldysh space. Sim- 
ilarly, A = (Ai,A2) T , and A stands for a 4 x 4 ma- 



trix A = [r+Ai - 
t± = (t x ± ir y )/2. 



r_A^](To + [t+A 2 - t-A^CTz, where 



In Eq. (2.1), Vo(q) is the bare static Coulomb potential 
between electrons in the dirty film. Below we will con- 
sider the situation when Vq (q) is screened due to the pres- 
ence of the nearby metal gate. In particular, if such a gate 
is situated at the distance b from the dirty film, parallel 

to it, V (q) = 27re 2 (l-e- 2fc 9)/ g ^5 47re 2 6. As it was dis- 
cussed at length in Ref. [GJ , the effect of long-range fluc- 
tuations of electromagnetic potential (which lead to the 
Coulomb zero-bias anomaly) is determined by the long- 
range part of the bare Coulomb interaction, and thus 
is suppressed once the Coulomb potential is screened, 
having the relative order of g^ 1 ln(87we 2 6) ln(l/f2r). We 
consider the case of a relatively short screening length 
b, and neglect long-range electric fluctuations and the 
Coulomb ZBA effects. On the other hand, Coulomb- 
induced repulsion in the Cooper channel Jl4|,[l5| does not 
depend on the long-range part of Vo(q) and is left un- 
changed by the screening gate. It is this effect of inter- 
action which we are going to consider in this paper. 



In terms of the a- model (2.1), diffusion-like collec- 
tive excitations of an electron system are described as 
slow fluctuations of the Q-matrix over the manifold 
Q 2 = 1. Small fluctuations near the metallic saddle 
point Q = A can be parametrized by the rotation ma- 
trix U = 1 + W/2 + . . ., with W obeying the relation 
{W, A} = 0. This constraint is resolved by 



W = u 



( W X T X + WyTy W + W Z T Z \ , . 

^ W +W Z T Z W x T x +WyT y J K ' 



where the matrix u is defined as 



1 F 

-1 



(2.4) 



K 



Variables with i = 0, z couple to diagonal in the 
Nambu space matrices and describe diffuson modes, 
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while Wj with j = x,y couple to off-diagonal in the 
Nambu space matrices and correspond to Cooperon 
modes. Their propagators are given by 

1 1 



(w i (q)e ie2 w l (-q) £2£l ) 

(wj(q) eie2 Wj(-q) e2£l ) 
(Wj(q) ei£a uJ,-(-q) e2ei ) 



iris Dq 2 
1 



1 



iris Dq 2 
1 



i[e\ 
1 



(2.5) 



7tz/ Dg 2 + i(ei + £2) 



General contraction rules for averaging of (Tr AW ■ 
Tr BW) over W are listed in App end ix |A|. 

In two dimensions, the action (2.1) can be studied by 
the renormalization group approach |fl5Ul| . At each step 
of the RG procedure one has to eliminate fast modes 
with either max(Z?g 2 , ei — e 2 ) > 0, (for diffusons) or 
max(Dg 2 ,e 1 + e 2 ) > f2 (for Cooperons) where 17 is the 
running ultraviolet RG cutoff. There are two types of 
logarithmic corrections to the action (2.1). Weak local- 
ization and interaction corrections have the relative or- 
der of g^ 1 ln(l/Or) and modify the conductance g of the 
system. Quantum corrections of the other type have the 
relative order of g~ x ln 2 (l/r2r) and emerge in renormal- 
ization of the Cooper channel coupling A. They become 
essential at such scales where localization corrections are 
still small and can be neglected. In this case A is the only 
running parameter of the bulk action which satisfies the 
RG equation |t!|§ 

-A 2 + A 2 „, A = ^, (2.6) 



2V5 



where £ = ln(l/f2r) is the logarithmic variable. The 
first term in the RHS of the RG Eq. (2.6) is due to the 
usual BCS logarithm, whereas the second contribution is 
the effect of Coulomb repulsion between slowly diffusing 
electrons. 



The bulk action (2.1) describes dynamics of the elec- 
tron system in the metal film. Possibility of tunneling 
between the island and the film is taken into account in 
the lowest tunneling Hamiltonian approximation by in- 
troducing the boundary term in the action: 



ITT J 



Tr r Qi s iQ. 



(2.7) 



Here Qi s \ refers to the island, spatial integration under 
the trace is taken over the area of the interface T, and 7 
is the dimensionless normal-state tunneling conductance 
per unit area of the boundary. 

Below we will consider thick enough superconductive 
island so that the absolute value |A| of the order param- 
eter is not suppressed by proximity to the normal film. 
The corresponding condition reads |A| 3> Gt/vVisI, with 
Vjisi being the island's volume. Also, the size of the island, 
d, is supposed to be smaller than the superconductive co- 
herence length so that for T <C T c the only relevant de- 
gree of freedom in the island is the phase <p of the order 



parameter. In the Keldysh formalism, one has to intro- 
duce its classical, ipi, and quantum, ip2, components, and 
arranging them into a matrix tp — <£>i<to + (p2&x one can 
write Qs in the subgap limit (e < A) as 



Qs = —ir + e lip + it- 



(2.8) 



We will also consider below the example of a normal (non- 
superconductive) island connected to the film by a tunnel 
barrier and biased at some voltage with respect to it. In 
that case the island's Q matrix can also be formally ex- 
pressed in terms of the superconductive phase: 



IN 



'■v/ 2 h_ e ~ iT *v/ 2 



(2.9) 



with ip/2 having the meaning of the single-particle phase 
which is conjugated to the number of electrons on the 
island (while ip is conjugated to the number of Cooper 
pairs). In both cases the applied voltage V(t) can be ac- 
counted by the Josephson relation dtpi/dt = 2eV for the 
classical component of the phase. 



III. FUNCTIONAL RG: NONINTERACTING 
CASE 



A. General idea 

In this section we will apply the renormalization group 
treatment to the system consisting of a small island cou- 
pled to a diffusive metal film by a tunnel barrier of ar- 
bitrary transparency. The effects of the Cooper inter- 
action in the film will be considered later on in Sec. 0, 
while now we will study noninteracting electrons. In this 
section we will restrict ourselves to the zero-energy limit 
when both temperature T, voltage eV and frequency 10 
are smaller that the Thouless energy = D/L 2 , and 
the perpendicular magnetic length Ih = y/^o/H ($0 is 
the flux quantum, H is the applied magnetic field) is 
shorter than the system size L. The effect of larger T, 
V or H is to destroy Cooperon coherence at scales larger 
than min{ y/D/T, y/D/eV, Ih}, and will be studied in 
Sec. VL The effect of high-frequency (ui ^> £rh) voltage 
is two-fold: it suppresses Cooperon coherence and also 
perturbs local electroneutrality of electron system. We 
will not study this effect in the present paper. 

The method developed in the present section can be 
applied both to the normal and superconductive island. 
However, since the rest of the paper will be devoted 
mainly to the study of charge transport in the N-S sys- 
tem, all formulas will be written for the case of a su- 
perconducting island, Q la i = Qs- In order to translate 
them for the case of a normal island, one has to substi- 
tute Qs by Qn, and the word "Cooperon" by the word 
"diffuson" . 
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FIG. 2. The second order in 7 correction to the action. The 
horizontal line schematically designates the N-S boundary F, 
with the two semicircles below it denoting Qs in the island, 
and the two-line object above it representing the Cooperon in 
the normal metal. 



Renormalization of the total action S'buik + Stun de- 
scribing an N-S contact is a complicated task. The main 
problem one encounters here is that this action does not 
reproduce itself under renormalization. In the lowest or- 
der over the tunneling transparency 7 studied in Ref . JO , 



the boundary term (2.7) generates the next-order term 
72 Ti(QsQ) 2 in the action, with 72 obeying the RG equa- 
tion <h)i I oc 7 2 . The physical meaning of this relation 
is simple. It describes the process when a Cooper pair 
tunnels from the island, coherently propagates in the nor- 
mal metal as a Cooperon and then returns back to the is- 
land. The resulting expression is bilinear in Qs (taken at 
the times of the first and second tunneling) and involves 
logarithmically (in 2D) large factor due to the probabil- 
ity of returning to the original point. The corresponding 
diagram is shown in Fig. H. 



b) 




agrams (a) and (b) conform well to this scheme, however 
it is not clear how to treat the bulk nonlinearity in the di- 
agram (c) . To overcome this problem we will employ our 
freedom to choose arbitrary parametrization of the field 
Q in terms of the matrix W, that can change the relative 
contribution of different diagrams. For the problem in 
question, the choice of the exponential parametrization, 
Q = Aexp(VF), nullifies all diagrams with Hikami boxes 
within the RG precision. To see this, we note that the 
leading logarithmic contribution arises from the integra- 
tion over fast momentum, so that the energy dependence 
of the diffusion propagators ( |2.5| ) can be neglected. To 
this accuracy, the matrices W and VW commute and 
the term Tr(VQ) 2 = — Tr(VVF) 2 becomes Gaussian in 
W , therefore there are no nonlinear vertices in the bulk 
in the zero-energy limit. Note also that the same expo- 
nential parametrization is used usually for the solution 
of the Usadel equations |25|,[l7j for the Green functions of 
disordered superconductors. 

Thus, in the exponential parametrization, diagrams 
of the type shown in Fig. |j|(c) are absent and one can 
carry out the RG procedure based on the tree-like struc- 
ture of the leading diagrams. Following this approach we 
present the total action as a sum of the bulk part S'buik 
and the part Sr which arises from the single tunneling 
action (|2.7l) under renormalization: 



Sr 



OO 

E 



S n 



2 ff ^7nTr r (Q s Qr. 



(3.1) 



c) 




FIG. 3. Three leading logarithmic corrections to the action 
at the fourth order in 7. 

In order to clear up the structure of the leading loga- 
rithmic corrections at arbitrary interface transparencies, 
let us consider the diagrams of the forth order over 7 
shown in Fig. ||. All of them ought to have the tree-like 
structure since any loop yields an additional small factor 
(/g -C 1. For two of them, labeled by (a) and (b), succes- 
sive Andreev tunnelings are connected by free Cooperons 
and diffusons, whereas the process shown in Fig. || (c) in- 
volves 4-Cooperon interaction vertex (Hikami box) in the 
bulk. Due to spatial integration, each of these diagrams 
evaluates to the third power £ 3 of the logarithm, and 
they exhaust the number of the most divergent diagrams 
of the forth order. 

Fig. U suggests an idea that the tree-like structure of 
corrections to the action may be suitable for RG ap- 
proach. Indeed, for each tree one can find the branch with 
the smallest momentum, say q. If we cut the whole tree 
over this branch, it will break up into two pieces which 
could be calculated at the previous steps of RG procedure 
since they contain only momenta greater that q. The di- 



The renormalized boundary action ( p.l[ ) describes pro- 
cesses of multiple Andreev tunneling coupled by coher- 
ent propagation of two-particle diffusion modes in the 
normal metal. It is determined by the infinite set of co- 
efficients ( "charges" ) j n and is a functional of the phase 
Tp of the island's order parameter. Once 7„'s are known, 
one can easily calculate the Andreev conductance and 
other quantities related to charge transfer through the 
system (see section [V B| ). Our main purpose below will 
be to establish the law of transformation of the coeffi- 
cients 7„ and to derive the corresponding multicharge 
renormalization group equation. 



B. Derivation of the multicharge RG equation 

According to the general approach discussed above, at 
each step of the RG procedure one has to calculate a cor- 
rection ASk, n to the action which results from c onn ecting 
two pieces Sk and S„ of the boundary action ( |3.l[) by a 
fast diffusion mode. This gives logarithmic contribution 
at scales larger than the island size d and, consequently, 
at energies SI smaller than u>d = D/d 2 . Therefore in this 
case it is convenient to define the logarithmic variable as 



C = ln 



(I 



2 In 



R 



(3.2) 
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where the current spatial scale R is related to Vl through 
il = D/R 2 . For the purposes of RG, the form of the 
island is of no importance, and all expressions depend 
only the total normal-state interface conductance Gt- 
Hence, one can formally consider a point-like island, with 



the transparency 7 in Eq. (2.7) becoming 7 = Gt and 
Tr r X = TrX(r = 0). 

Our aim is to derive an effective low-energy action of 
the phase Tp on the island, therefore we consider Qs as 
a slow variable and decompose only metallic Q into the 
fast and slow parts. Such a decomposition consistent 
with the exponential parametrization adopted is given by 
Q = U~ 1 Aexp(W')U, where U is a slow rotation matrix 
and W' is a fast variable Jl],|l5| . The relevant interaction 
vertex reads 



-nr 2 gn ln Tr r (Q s Q) n - 1 QsU- 1 AW'U, (3.3) 



where the additional factor of n accounts the fact that 
the fast diffusion mode can be coupled to any matrix Q in 
S n , and Q = U~ 1 AU. The form of the vertex cor- 
responds to the tree-like structure of the diagrammatic 
expansion discussed in Sec. Ill A. The lack of higher- 
order terms in W in Si n t,n amounts to neglecting closed 
loops in virtue of the small parameter Q/ g. 

For k 7^ n, the correction to the action is given by 
A5fe,„ = i(Sint,fc£int,n)- This average (in the zero-energy 
limit) is given by Eq. ( |A3[ ), with 

(3.4a) 
(3.4b) 



A = U(Q s Q) k ~ 1 Q s U~ 1 A, 
B = U{Q s Q) n - l QsU- l k. 



Calculating the logarithmic integral over the fast momen- 
tum and omitting the tilde sign over slow variables one 
obtains 

AS*,-,™ = ^y^fcn7 fc 7„AC 
x Tr r [(QsQ) k+n ~ {Q sQ)^ 1 Q s{Q sQT' 1 Q s] ■ (3.5) 

Employing invariance of the trace under cyclic permuta- 
tions and 
for k =/= n 

™ 9,. v tv r„, /.,/-» _ (Q 5 Q)l*-«lj. 

(3.6) 



tions and using the relations Q 2 = Q 2 S = 1, we get finally 



ASk.n — 



- kn lkln AC Tr r [(QsQ) k 



Thus we see that averaging of the terms Sk and S n 
from the action (3.1) modifies the terms Sk+n arL d 
S\k- n \ witn A7 fe+ „ = -(l/2)fcn7fc7„AC and A-y| fc _ Tl | = 
(l/2)fcn 7fc7n AC. 

Analogous calculation of ASk,k — (i/2)(S? nt k ) yields 

A5 fc , fc = ^A ; 2 7, 2 ACTr r (g s g) 2fc (3.7) 

that results in renormalization of the coefficient A^k = 
-(l/4)fc 2 7, 2 AC. 



Taking into account all possible pairings ASk, n be- 
tween the terms of the action (3.1), we derive the main 
equation that describes evolution of the charges j n under 
renormalization: 



\ ^2 ^ ~ fc ) ^ln~k + \^2K n + k ) lkln+k- 

(3.8) 



k=l 



fe=l 



Initial conditions for these multicharge RG equations at 
C = follow from the bare tunneling action (2.7): 



7n>l(0) =0. 



(3.9) 



Note that Q 2 S — 1 is the only property of the is- 
land's Qs which was used to derive the RG equations 
( j3.8| ). Therefore these equations universally describe 
both normal-metal and superconductive islands on the 
same footing. 

Equations (3.£) look quite complicated. Their hid- 



den algebraic structure becomes transparent after Fourier 
transformation. To this end we introduce a 27r-periodic 
function u(x) of an auxiliary continuous variable x ac- 
cording to the definition 



u{x) 



n^n smnx. 



(3.10) 



Transforming Eq. (3.8) into x-representation we obtain 
the following RG equation for the function u(x): 



= 0. 



(3.11) 



The differential equation (3.11) is nothing but the well 
known Euler equation describing ID motion of a com- 
pressible gas, with £ playing the role of time. The 
initial condition for Eq. ( |3.11 ) at £ = is given by 
uq(x) = a sinx. 



The Euler equation (3.11) can be easily solved with the 
help of characteristics. The latter are described by the 
following system of differential equations: 



dx 



du 



0. 



(3.12) 



They are trivially integrated: u — uq(xo) and x = xq + 
uq(xq)(. In order to get u(x) = uq(xq(x)) one has just 
to find the inverse function xq(x). For uq(x) — a sin a;, 
the function u(x) is implicitly defined by the relation 



u(x) — asin(a; — u(x)(). 



(3.13) 
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FIG. 4. Solutions for the Euler equation ( 3.11 ) at different 
values of t = a£. 

The solution of the Euler equation is conveniently ex- 
pressed through the variable t — a( — Gt/Gd which is 
equal to the ratio of the tunneling conductance of the bar- 
rier Gt to the conductance of the metal film Go — 4ng/C 
at the current scale. The solution for different values of 
t is shown in Fig. The derivative u x (ir) diverges at 
t = 1, exactly at the point where Gd = Gt- For t > 1, 
the discontinuity ("shock wave") develops at x = ±tt. 
For t ^> 1, the solution acquires a saw-like type with 
u(x) ~ (a/t)x = (Go/4:^g)x. 



the island and in the lead. Below we will assume that 
the metal lead is at zero bias, so that Qn reduces to A 



and the Proximity Action (3.14) acquires the form 



-'prox 



oo 



7„Tr(Q s A) f 



(3.15) 



In this formula, 7„ are taken at "time" t = aCTh, with 
(/rh = \n{u)d/ETh) = 2\n(L/d). It can be calculated 
through the integration of the FRG equation ( |3.11 ) for 
the function u{x, (,) up to the scale (Th- 
in order to express the action in terms of the phase 
~ip, one has to substitute Qs from Eq. (2.8). Taking the 
trace over the Nambu space, we obtain that only terms 
with even n contribute to the Proximity Action: 

oo 

S pm *lQs} = -2i7r 2 5 ^(-l)"72„Tr K (e^Ao e -^A )", 

71=1 

(3.16) 

where the subscript 'K' indicates that the matrix trace 
is t aken o nly over t he Ke ldysh space. As explained in 
Sec. IIIB , the action ( |3.15 ), after substitution Qs — > Qn, 
also applies to the case of the normal island. Using 
Eq. (2.9) for Qn and taking the trace over the Nambu 
space, we arrive at the following expression for the action: 



C. Proximity Action 

In the previous section we have established the trans- 
f or mat i on l aw for the coefficients 7„ of the multicharge 
action (3.1) under the action of the renormalization 
group. After eliminating diffusion modes up to some 
value of the logarithmic variable the action becomes 
a functional of the island's Qs and Q{R) in the film 
taken at the special scale R ~ de^ 2 from the island. 
For a finite system bounded by a perfect metallic lead 
located at a distance L (L ^> d) from the island (see 
Fig. 0) , renormalization of the action started at the en- 
ergy scale cod should stop at the scale of the Thouless 
energy £rh = D/L 2 when all diffusive electronic degrees 
of freedom in the dirty metal are integrated out. The 
resulting action, 



■SproxIQs, Qn] = 



-iTT 2 g 



n=l 



Ty(Q s QnT, (3.14) 



thus becomes the functional of the island's Qs, Eq. (2 



and the matrix Qn, Eq. (2.£), in the lead. Such an action 
resulting from elimination of the modes in the diffusive 
conductor will be referred to as the Proximity Action. 

According to Eqs. (2.8) and (2.£), the only degree of 
freedom in the matrices Qs and Qn is the one associ- 
ated with the phase rotation. In other words, externally 
controlled voltage is the only dynamical variable both in 



S[Q N ] = -2iir 2 g V 7n TrK(e^/ 2 A e-^/ 2 A )™. (3.17) 



n=l 



Remarkably, our function u{x, Ctii) can be directly re- 
lated to the generating function of transmission coeffi- 
cients F(<f>) introduced by Nazarov [Q (cf. also Ref. p2|). 
Comparing our Eq. ( |3.13| ) with Nazarov's Eq. (19) one 
can verify that 



F{4>) 



Ang u((f>) 



sin < 



(3.18) 



Note also that appearance of discontinuity for t > 1 
goes in parallel with opening of conducting channels with 
transparencies T — > 1 Q. The saw- like solution for the 
function u(x) at t 3> 1 results j| in the Dorokhov's |2(| 
bimodal distribution function V(T) = G D /2T^T^T for 
the transmission coefficients of a diffusive system. 

It might be surprising to realize equivalence of the two 
approaches since Nazarov's result is valid for any ge- 
ometry while our FRG equation apparently holds only 
in 2D. Such a coincidence is related to the fact that 
both Nazarov's derivation and our derivation of the func- 



tional RG equation (3.11) implied the zero- energy limit, 
T, eV, us < E^h ■ Under this condition it is possible to 
derive the universal distribution of the eigenvalues of the 
transmission matrix as well as to neglect Cooperon de- 
coherence in deriving Eq. (3.11). Note further that sum- 
mation of the tree-like diagrams for the action can be 
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performed without the use of the RG at all: one has to 
integrate diffusion propagators over momentum indepen- 
dently for each branch of the tree. But in the zero-energy 
limit, such integration is expressed in terms of the Green 
function of the Laplace operator, thus reducing to the 
total resistance of the system. Therefore, in this case 
Eq. ( |3.11 ) rewritten in terms of t = Gt/Gd remains 
valid for any geometry of the disordered normal region. 

Power of the Keldysh multicharge Proximity Action 
will be crucial if one has to go beyond the zero-energy 
limit and, especially, for the problems with interaction. 
The corresponding FRG equ ations replacing the zero- 
energy-limit Eulcr equation ( |3.1l| ) will be derived in 
Sees. and 



IV. PHYSICAL QUANTITIES 



Now we will discuss how the action ( |3.15| ) together 
with the FRG equation (3.11) valid in the zero-energy 
limit can be used to calculate physically observable quan- 
tities. In this manuscript we will focus on the conducting 
properties of the system between the island and the lead. 

Once the Proximity Action S prox [Qs] is known, the 
system conductance as well as current statistics can be 
calculated following the method of Ref. 0] . We will sup- 
pose that the island is connected to the voltage source 
by an ideal conductor so that phase fluctuations are sup- 
pressed (the effect of phase fluctuations will be stud- 
ied elsewhere [P7| ). In this case the general expres- 
sion Z = J D(piDif2 exp{iS p[ox [ipi,ip2]} for the Keldysh 
partition function reduces to Z = exp{iS pmx [ip s i, ip S 2]} 
where ~ip s is the island's phase controlled by the source. 
Dynamics of its classical component is governed by the 
applied voltage through the Josephson relation: 



dips 



df 



2eV(t), 



(4.1) 



while the current operator is given by the derivative with 
respect to the quantum component of the phase: 



lit) 



(4.2) 



In order to obtain the expectation value of currents taken 
at different moments of time one should apply the oper- 



ators (4.2) to InZ: 



(1(h) . . . I(t k )) = 1(h) . . . I(t k ) lnZ[<p el , <psi] 



fs2=0' 

(4.3) 



Below we will employ these relations to calculate the 
conductance and noise for the case of the normal and 
superconductive islands. 



A. Normal junction: conductance and noise 

To illustrate the formalism, we start by considering the 
conductance of the system for the case of the normal is- 
land. We will trace how addition of resistances Rt and 
Rd is realized within the multicharge action approach 
and calculate the current noise power. 

The current response (I(t)) = Y^=i(I(t))n to the bi- 
ased voltage V(t) is given by the sum of the partial con- 
tributions arising from the terms of the action (3.17): 



(*(*)>» 



o . 2 6Tr K (e^/ 2 A e-^/ 2 A a y 
2ien g j n 



5<p 2 (t) 



ip 2 =0 



(4.4) 



(here and below we om it the source phase subscript 
's'). Now using Eq. (B5) we calculate the current / = 
(e 2 /h)GV, where the dimensionless conductance G is 
given by 



G 



n=l 



(4.5) 



In terms of the function u(x) it reads 
G = 47rgu x (0). 



(4.6) 



According to the Euler equation ( 3.11 ), G obeys 
d(l/G) 1 



(4.7) 



d{ 4tt 5 ' 

that amounts to the anticipated addition of resistances 
Gt GtGd 



G 



1 + i Gi 



G, 



where 



Gd 



2ng 
\n(L/d) 



(4.8) 



(4.9) 



is the metal film conductance between the island and the 
lead. 



Eq. (4.8), though derived in the low-frequency limit, 
u) <C D/L 2 , can be shown to remain valid for much larger 
frequencies, up to w m ax = 2na/L (for the bare Coulomb 
interaction) or w max = Aita/b (if the Coulomb interac- 
tion is screened by the presence of a gate situated at the 
distance b from the film, see Sec. [n]). This frequency 
determines the time scale w~^x a ^ which the electroneu- 
trality of the system sets in. 

After this lesson let us turn to a more involved calcula- 
tion of the current-current correlator at a fixed constant 
bias V . Now one has to apply two current operators 
to the action (3.17). Similarly to the above consider- 



ation we present the Fourier transformed correlator as 



UJ- 



with 
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2e 2 7r 2 3 7, 



8<p 2 (w) 5(f2(—U)) 



V>2=0 

(4.10) 



The corresponding functional derivative is calculated in 
Appendi x p3| . Substituting it from Eq. (B15) and employ- 
ing Eq. ( E-q ) one obtains 

(J w J- w ) = ^{(3-P Ar (t))*H 

+ \PN{t) [*(w - eVO + #(w + eF)]}, (4.11) 
with the function \&(u;) defined as 

poo 

= y [l-F(£+)F(£_)]d£ = u;coth— , (4.12) 
where E± — E ±lj/2. The function Pjv(t) is given by 

2Wa:a;^ (0) 



Pjv(t) 



1+&^ = 1 



u x (0) 



(4.13) 



The third derivative, ti xx:c (0), can be easily calcu- 
lated with the help of characteristics. The one reaching 
the point x = is obviously x(() — with xo = 0. 
Writing u{x) = uo(xo) and using the relation d x = 
[1 + t cos Xq] d Xo at Xq = we obtain 



and, therefore, 



(i+ty 



1 G 4 

3 
'J' 



Airg G? 



P N (t) = 1 + 2^. 

KJTrp 



(4.14) 



(4.15) 



Eqs. (4.11) and (4.15) provide a general description 
of the noise power for arbitrary relation between the 
frequency (which should be smaller than i?Th)j volt- 
age and temperature, and arbitrary Rt/Rd- Though 
it was derived in the zero-energy limit, the condition 
(T, eV) < -&rh can be relaxed since for a normal island, 
charge propagation in the diffusive conductor is described 
in terms of diffusons which are insensitive to decoherence. 



In the limit V — > 0, Eq. ( 4.11 ) reduces to the Nyquist- 
Johnson equilibrium thermal noise: 



it T \ z 2 G , w 



(4.16) 



In the limit eV ^> (lu,T), one obtains for the shot noise 
power: 



(4.17) 



that coincides with Nazarov's result Q obtained via the 
distribution function of transmission coefficients. The 



function P/v(t) is plotted in Fig. ||by the dashed line. For 
f > 1, when the system resistance is dominated by the 
diffusive conductor, the shot noise power is three times 
smaller than its Poisson value. This result was first ob- 
tained in Ref. [Q with the help of the Dorokhov's bi- 
modal distribution J2(|, and in Ref. j29), in the frame- 
work of the classical Boltzmann equation with Langevin 
sources. 



The Pr oxim ity Action ( 3.15 ) together with the FRG 
equation (3.11) allows, in principle, for the calculation of 
higher momenta of current fluctuations and even the full 
statistics of transmitted charge (cf. Refs. |3(],|3l|] ) . We 
leave this problem for future investigation. 



B. Andreev conductance 

In this section we will turn to the case of a supercon- 
ductive island. We will assume that the island is not very 
small so that its order parameter satisfies A > uj^. This 
means that at the scales < Ud relevant for RG, quasi- 
particle reflection from the island is of A ndreev type and 
one can use the subgap expression (|2.8|) for the island's 

Qs- 

Similarly to the N-island case, we write the current as 
a sum of the partial contributions (I(t)) n from different 



terms of the action (3.16) 



</(*)>„ = {-l) n 2ie7T 2 g l2n 



l 5Tr i f(e^A e-^Ao) r 



Stpz(t) 



Calculating the derivative with the help of Eq. 
find for the dimcnsionlcss Andreev conductance: 



(4.18) 



) we 



G^ = 16tt 5 ^(-1) 



n 72n, 



or, in terms of u(x): 



G A = ^gu x {\). 



(4.19) 



(4.20) 



Note that Eqs. ( [1.6] ) and ( i.2C| ) for the case of the normal 
and superconducting island look very similar. The only 
difference is at the point where the derivative should be 
taken. 

Now we will use the general expression ( |4.2CD to 
calculate G a in the noninteracting case described by 



Eq. (3.11). The simplest way to calculate u x {^) is to 



use characteristics of the Euler equation. To do this, one 
has to find the one which leads to x — ir/2 at "time" t. 
Writing its initial point as xq = tt/2 — 0(t) we obtain the 
following equation for the function Q(t): 



6(t) =*cos6(t). 



(4.21) 
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For small t, 0(f < 1) = t + 0(t 3 ) while in the limit of 
large t, 0(f > 1) = tt/2 - n/2t + 0(f~ 2 ). Now, applying 
d x = [1 + tcosxo] d Xo to it(x) = uq(xq) we obtain 



*s( 2 ) 



a sin 0(f) 
1 + f sin 0(f) 



and, finally, 



G 



G 7 



sin 0(f) 
1 + f sin 0(f) 



= G 



D 



f sin 0(f) 
1 + f sin 0(f)' 



(4.22) 



(4.23) 



where Gd is defined in Eq. (4.9) 



In terms of the resistance, the zero-energy result (4.23) 
can be written as 



Ra = Rd + Rr,eS, 



(4.24) 



where i?T,cff(f) = Rt/ sin 0(f) has the meaning of an ef- 
fective interface resistance connected in series with the 
metal-film resistance Rd = (h/e^G^ 1 . The former de- 
creases with the increase of Rd due to disorder- induced 
enhancement of Andreev reflection determined by the 
probability of ret urning to the origin. 

Equation ( 4.23 ) for the Andreev conductance had been 
previously derived by Nazarov both with the help of 
the Usadel equation || and the distribution function 
F(<f>) of transmission eigenvalues 0. His latter approach, 
though formulated in a very different language, appears 
to be directly connected to our multicharge action ap- 
proach. Using the relation between F(<j>) and u(x) es- 
tablished in Eq. ( |3.18] ), one can verify that Nazarov's re- 
sult |] G a = (<9i ;l /90)|0 =7r /2 exactly coincides with our 
Eq. (120|). 



C. Noise of NS current 



Decomposing the current-current correlator into the 
sum of the partial contributions, (I u I-J) n , similarly to 
Eqs. (4.18 ) and ( 4.10 ), using Eq. QB14| ), and employing 
Eq. (4.19) we obtain 



e 2 G A 
3h 



{(3-P s (t))*(w) 



^Ps(f) - 2eV) + V(lu + 2eV)]j, (4.25) 



where ^(w) is defined in Eq. (4.12). The superconductive 
noise function Ps(f) is given by 

P^) = l + 2 ^ ( "5f 72 " =l-^#. (4-26) 



E„=i(-i)"« 2 7 2 , 



2« a (f) 



Eq. (4.25) has the same form as its analogue for the N 
island with the replacement P/v(t) — > Ps(f) and eV — > 
2eV. 



N-S 



N V 



_l I I L_ 



_l I I L_ 



_l I I l_ 



0.0 



0.5 



t 



1.0 



1.5 



FIG. 5. Noise functions Pjv(f) (dashed line) and Fs(f) 
(solid line) for normal and N-S junctions, correspondingly, 
vs. the ratio t = Gt/Gd- 

For the noninteracting system, the value of u x (^) was 
calculated in Eq. ( 4.22 ), and calculating u xxx (^) in the 
same way we obtain 



U XXX ( o ) 



Thus, we get 



a sin 0(f) 



3af cos 2 0(f) 



(1 + fsin0(f)) 4 (1 + fsin0(f)) 5 ' 



(4.27) 



Ps(t) = 1 



1 + tan + 30 cot 



2(1 + tan©) 4 



(4.28) 



where = 0(f) is defined in Eq. (4.21 



In the li mit o f vanishing bias, Eq. (4.25) with Ps{t) 
from Eq. ( 4.28 ) reprod uces the Nyquist- Johnson ther- 
mal noise giv en by Eq. ( 4.16 ), with G being substituted 
by Ga, Eq. ( 4.23 ). In the opposite limit of large bias, 
eV ^> (oj,T), it gives the power of the shot noise be- 
tween the normal and superconducting terminals ([Tsf : 



(II) 



shot 



2e 2 G A 
3h 



Ps(t) eV. 



(4.29) 



The function Ps(f) is plotted in Fig. || together with its 
analogue for the normal junction Pjy(f). On increasing 
f, both Pjv(f)/3 and Ps(f)/3 evolve from the value of 1 
describing the Poisson noise at a tunnel contact to 1/3 
corresponding to the diffusive conductor. An excess fac- 
tor of 2 in Eq. ( 4.29| ) accounts the fact that in an N-S 
system charge is transferred by Cooper pairs. 

and IV C| refer to the 



Results obtained in Sees 
zero-energy limit 



[VB 



In Sec. VI we will extend our theory 



to the case of large temperature or voltage, max(T, eV) > 
Em- 
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V. EFFECT OF INTERACTION IN THE 
COOPER CHANNEL 

A. General idea 

Here we calculate the effect of repulsion in the Cooper 
channel on the renormalization of the effective boundary 
action (3.1). In the lowest order over the interface trans- 



parency, the RG equation for 71 was derived in Ref. 
We will extend this result to account the whole set of 
charges {7„}. 



As in Sec. [II, we first turn to the analysis of the lowest 



order perturbative corrections to the action. The Cooper 
channel interaction is described by the following vertex 
in the bulk action: 



S\ — 



J dvdttxa x [Q 2 - (t z Q) 2 



(5.1) 



obtained from the a- model action (2.1) by eliminating 
the fluctuating A field [Q . The tree-like diagrams for the 
action of noninteracting system should be modified now 
by all possible insertions of the vertexes S\. In the low- 
est order in the interface transparency 7 (cf . Fig. |^) , that 
was carried out in Ref. Q|. The diagrams of the forth 
order in 7 and the first order in A are shown in Fig. |^. 
They are obtained from the diagrams of Fig. || by inser- 
tions of the Cooper vertex denoted by a dot. The latter 
cuts a diagram into two parts with independent energy 
integration variables. 





FIG. 6. Leading diagrams in the forth order in 7 and the 
first order in A. Shadowed is the piece of the diagram where 
the logarithmic contribution of the Cooperon adjacent to the 
A- vertex results from the energy integration. 

The leading logarithmic contribution of a diagram from 
Fig. ^ is due to the momentum integration in each dif- 
fusion mode except for one of the Cooperons adjacent 
to the vertex S\, whose contribution is logarithmic due 
to the energy integral J dE/E. The corresponding part 
of the diagram containing such a Cooperon is shadowed 
in Fig. ||. Thus, each of the diagrams shown in Fig. ^ 
evaluates to the forth power of the logarithm, cx \j 4 ( 4 . 



To achieve this situation, the energy E and momenta 
in the shadowed part of the diagram must obey some 
inequalities. First, the momentum of the Cooperon ad- 
jacent to the A-vertex, q, must satisfy the relation Dq 2 < 
E. Second, momenta of all other shadowed Cooperons, 
qi, must satisfy the relation Dq 2 > E. Only under these 
conditions the diagram evaluates to the maximal power 
of the logarithmic variable C- 

Were all diffusion modes in the shadowed region of the 
diagram Cooperons, the energy E would be the slow- 
est variable in this region. Then, on the level of the 
RG, the shadowed region would be taken into account 
by the pairing of the vertex S\ with the multicharge ac- 
tion Sr [Eq. ( |3.lD ] obtained at the previous steps of the 
RG procedure. Since E would be the slowest variable 
in the shadowed part, all metallic Q's entering the ac- 
tion Sr paired to S\ should be set to A indicating that 
no further pairings with slower variables are allowed to 
the shadowed region. As a result, as will be shown in 



Sec. VB, the shadowed region will give a correction to 



the term Si cx 71 Trr QsQ- However, the shadowed re- 
gion might contain diffusons as well. Their momenta qi 
are not restricted by the inequality Dq 2 > E and may 
be "slower" than the energy E. Were this the case, the 
multicharge action ( |j.l| ) would not reproduce itself upon 
pairing with S\ since one would not be able to perform 
the integral over E as the result would depend on next 
steps of the RG. 

Fortunately, diffuson contribution in the shadowed 
part of a diagram can be disregarded. This can be seen 
already from the diagrams of the fourth order shown 
in Fig. |(| Among them only the diagram (c) con- 
tains a shadowed diffuson. However, direct calculation 
shows that its contribution to the action vanishes identi- 
cally [p2| . This statement can be generalized to arbitrary 
order. One can prove that any diagram containing a dif- 
fuson in the shadowed region is exactly zero. Physically, 
the shadowed part of a diagram accounts the correction 
to the Usadel spectral angle 8 (cf. Appendix ^|) which is 
determined solely by Cooperon modes. 



B. Functional RG equation 

To calculate the correction to the action due to the 
pairing S n ,\ — i(Si n t,nSint,\) one should extract one fast 
W from S x : 



ti 2 v\ 



5int,A 



Tr{a x ,Q}A(W' -t z W't z ), (5.2) 



and another one from S n according to Eq. (3.3). The 
only fast variable to be integrated out at this step of the 
RG procedure is the energy E running over the Cooperon 
propagator (WW). This fast energy runs also over all 
Q's under the trace in Smt.n- Substituting Q = A(E) 
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into Eq. (f37!|) as explained above and performing the av- 
eraging ove r th e fast Cooperon in {S- m t tn S- m t t \) with the 
use of Eq. ( |A2D , we obtain that only the terms with odd 
n give a nonzero contribution to the action (hereafter the 
tilde sign over slow variables is omitted): 



S n ,\ = 



dr 



dt dE' dE 



xtv{(Q s A (E)Y 



i .; (2tt) 2 E 
(A (E)Qs) n ,a x }Q. 



(5.3) 



Since E is a fast variable, all Qs(t) are taken at the same 
time t, and Q = J" e tu]t Q E , _ u /2tE , +u) /2 (duj /2tt). The sign 
in the above equation comes from commutation of t z con- 
tained in A = Aqt z with Qs in order to transform (QgA) n 
to (Q s A ) n . 
Let us calculate 



M n = {(QsA ) n + (A Q s ) n ,a x }, 



(5.4) 



which enters Eq. ( |5.3| ). It can be written as .M, 
X~ n , <r x }, where X = Q S A . According to Eq. (|2.8|), Qs 
contains the classical and quantum parts: Qs = P + q&x- 
Then it is easy to show that M\ = 4QsF(E). Now 
we will show by induction that for all odd n, A4„ = A4\. 
Consider the difference M. n+ 2 — M. n - It contains X n+2 — 

X n _ x -n + X -n-2 = (jn+1 _ X^ 1 )^ - AT" 1 ) = 

T{X)(X - X- 1 ) 2 . Within the RG accuracy the distri- 
bution function F(E) = sgnE so that F 2 = 1. Under 
this condition, (X - A -1 ) 2 = [Qs, A ] 2 = that proves 
the statement M n = Mi- 
Substituting 

{(QsA (E)) n + (A {E)Q s ) n ,a x } = 4Q s sgnE (5.5) 



into Eq. (5.3), we get 

S n ,x = (-l) {n - 1)/2 iir 3 gnj n XA{TrrQsQ- 



(5.6) 



Therefore, only Si is subject to renormalization by the 
A-term: 



/ due to A 



-A]T(-1) ; (2; + 1)7 2J+1 . (5.7) 



1=0 



As a result, an additi onal t erm will appear in the RHS 
of the Euler equation ( |3 .11 ): 



-X(C) it(§ ) sin a;. 



(5.8) 



Here A(C) is described by the bulk RG equation ( |2.6| ) 
(tunnel coupling between the superconductive island and 
the film, though induces finite A in the bulk, does not 



dSx 
d( 



dSu 



have logarithmic effects on the Cooper interaction). Rs 
solution |l5|Jl|] is given by 



A(C) = 



\d + X g tanh A S C 

X~d ' 

1 + — tanh X g ( 



(5.9) 



with X e i defined at the energy scale u>d (cf. the definition 
.2|) of the logarithmic variable Q) : 



X n + A g tanh I X g In 



1 + -r^ tanh ( X In ■ 



A„ 



(5.10) 



and X n is the interaction constant at the energy scale 
H/t. 



C. Andreev conductance 

Here we consider the effect of repulsive interaction in 
the Cooper channel of a normal conductor on the An- 
dreev conductance, in two limiting cases of weak and 
strong interaction. Of course, we always assume that di- 
mensionless coupling constant A <C 1. We will see below 
that the limit of strong interaction means A 3> 
and thus can be realized at relatively poor transparency 
of the interface. However, we will see also that even the 
weak-A solution produces the interaction-induced correc- 
tion to Ga which grows logarithmically with the spatial 
scale L and eventually crosses over to the strong-coupling 
limit. 



1. First order correction 

In the limit of small A(£), it can be taken into account 
perturbatively. The equations for the characteristics of 
Eq. (5.8) have the form: 



dx 
d( 



du 

— =-A(C)u(f,C)sinz. 



(5.11) 



The presence of u(f , C) m Eq- (5-11) which is determined 
by the trajectory reaching x = 7r/2 in "time" Q intro- 
duces a nonlocal in £ coupling between different trajec- 
tories. In the lowest order over A(£), we search for the 
solution in the form x{C) — xq + a^sincco + &c(C) an d 
u(£) = asinoio + 5u(C)i where Sx and Su are propor- 
tional to A and vanish at £ = 0. Substituting this into 



Eqs. (5.11) and keeping only terms linear in A, we obtain 

(5.12) 



—— = — A(C) acos[0(a£)] sin[a;o + aCsin(xo)] 
etc, 



These equations can be easily integrated as 
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5x(Q = —a I — rj)\(rj) cos[Q (arj)] sin[xo + ar] sin xo]drj, 



HO 



-a / X(r]) cos[Q(ar])]sin[xo + ar]smxo]dr]. 
Jo 



(5.13) 
(5.14) 



Thereby, we know the characteristics of Eq. ( p.q ) in the first order over A(£). As we have already seen in Sec. [VB, 
in order to calculate the Andreev conductance Ga — 4^gu x (^), one has to search for the characteristic which 
leads to x — ir/2 in "time" £. Resolving equation x(Q — 7r/2 we find that the initial point of such a trajectory, 
x = 7r/2 — O(aC) + gets shifted by 



1 + a(sinO(a() Jo 



(C — r])\(r]) cos[6(a?7)] cos 



e«) i 



di) 



(5.15) 



compared to the noninteracting case. A nonzero A influences the derivative u K (§) = (du/dxo) / (dx/ dxo) which 
determines the Andreev conductance in two ways. The first correction, S^u x (j), to the quasiclassical result (4.22) 
is related to the change of the initial point xq of the trajectory. The second c orrec tion, 5^- 2 'u x {%), is due to the 
modification of the functional dependences of x(() and u(() versus xq by the terms ( [5. 13 ) and ( 5.14 ). A straightforward 
calculation yields 



and 



[l + <sine(aC)] 2 J 



S W U (*) - acosQ(aQ 

° Ux[2 >- [l + aCsinOK)] 2 



X(rj) cos[0(a?7)] cos 



e(ac) [ i - j 



[1 + a77sin9(aC)] 2 ^. 



(5.16) 



(5.17) 



Finally, with the help of Eq. ( 4.20 ) we obtain for the Andreev conductance in the lowest order over the Cooper- 
channel interaction A(£): 



Q(t) 



G A tsinO(t) 
G~d ~ 1 + <sin9(t) ~ [1 + tsin6(i)] ; 



■ C /" X(xO (1 - x) cos[0(t)(l - x)] dx 



:\. [1 + xtsinQjt)} 2 9(xt) . . ... 

C / A(afC) n , 77: r>^M2 sm[9(t)(l - x)} dx. 



ii 



[1 + tsine{t)} 2 x 



(5.18) 



It is interesting to compare Eq . (|5 .18 ) with the result 
of the standard (cf., e. g., Ref. |l7|) calculation of Ga 
by means of direct solution of the Usade l equ ations , pre - 
sented in Appendix]^. While both Eqs. (|5.18|) and (C22) 
contain corrections linear in A, their meaning is some- 



what different. Whereas Eq. ( 5.18 ) contains, in general, 
the renormalized coupling constant A(£) which depends 
on the energy scale through Eq. (5.9), the result ( C22 ) is 
expressed in terms of the bare coupling constant A^ de- 
fined at the scale u>d, cf. Eq. ( 5.1C ). Note that the Usadel 
equation contains, in principle, logarithmic dependence 
of A(C) which follows from the first term of Eq. (2.6): 
A(C) = Ad/(1 + XdQ- It emerges as a result of the self- 
consistent determination of the proximity-induced A in 
the normal metal, cf. Eq. (C3). However, in order to take 
these effects into account within the approach presented 
in Appendix it would be necessary to solve the whole 
problem with higher-order accuracy in A^, that seems to 
be an extremely complicated task. The main compli- 



cation is the appearance of the supercurrent flowing in 
the normal region that leads to the rotation of the phase 
<p(r) of the anomalous Green function and also mixes 
both components, / and /i, of the distribution function. 
As a consequence, the problem becomes hardly tractable. 
Staying within first-order accuracy over A<j, we identify 
A(C) = Ad and find exact agreement between Eqs. ( |5.18| ) 



and (C22). Moreover, two c ontrib ution s to the correction 
to Ga given by Eqs. ( 5.16 ) and ( 5.17 ) have their direct 
counterparts within the Usadel method: the first one is 
due to th e modification of the spectral angle 6(r = d) [cf. 
Eq. ( |C10] )], whereas the second one is, physically, due to 
the correction to the distribution function fi(r = L) by 
the presence of the induced A [cf. Eq. (C21)]. 



Another principal limitation of the quasiclassical ap- 
proach based on the Usadel equation is that the latter 
does not allow for an account of fluctuation corrections, 
such as, e. g., Finkel'stein correction to A(C). This is re- 
lated to the fact that the Usadel equations are nothing 
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more than the mean field approximation for the Keldysh 
action ( |2.l| ). Therefore, our approach seems to be the 
only tool for studing the effect of the Cooper channel 
interaction on the conducting properties of 2D N-S sys- 
tems. 



A somewhat cumbersome expression (5. IS) can be 
simplified in the limiting cases of small and large t. 
For t -C 1, the noninteracting Andreev conductance 



G 



(0) 



Gji/Gd, and one has for the interaction correc- 



tion 5Ga = Ga — G^a' 



SG a 
G a 



-2C / A«)(l -x)dx. 
'o 



(5.19a) 



In the opposite limit, t > 1, the conductance is deter- 
mined by the diffusive conductor, G? ~ Go, and 



SG A 
G a 



w 7TX 

X(xQ x cos — ax. 
o 2 



(5.19b) 



An important consequence of Eqs. (5.19) is that the rela- 
tive value of the A-correction to the subgap conductance 
is proportional to AC at any value of Gt/Gd, includ- 
ing the case of a completely transparent interface with 
G" 1 = 0. Thus, our result differs from the one ob- 
tained by Stoof and Nazarov jl7|] who considered the 
same kind of problem in the ID geometry and found a 
weak regular correction ~ A only. The growth of the 
interaction-induced correction in 2D indicates that qual- 
itative changes in the behavior of the function u(x) do 
occur at large enough scales C and one has to sum the 
leading logarithmic terms, (AC)™. 

In the most general case of arbitrary A(C), the FRG 
Eq. (5.8) should be solved numerically. However, in the 
limiting cases of strong and weak Cooper interaction, an 
analytical solution can be obtained. The forme r refers to 
the situation when the first correction, ( 5.19a ), becomes 
of the order of one already in the limit of weak trans- 
parency of the barrier, ( « 1. The latter applies when 
the first correction, ( [5. 19b ), becomes important only in 
the diffusive regime, f > 1. The boundary between the 
limits of strong and weak repulsion in the Cooper chan- 
nel slightly depends on the relation between Ad and X g . 
In the rest of this section we will consider the case when 
Ad ~ X g . Then repulsion is strong (weak) provided that 
X g 3> a (X g -C a). We want to emphasize that the no- 
tion of strong/weak interaction refer not to the value of 
A itself but to the ratio A/a = AirgX/Gx- In both cases 
the interaction correction can be either small or large 
depending on the value of AC = ingX/Go- 

The limiting cases of strong and weak repulsion will be 
considered Sees. V C 2 and V C 3, and then (in Sec. V C 4) 
we will turn to the discussion of the general case of arbi- 
trary A. 



Eqs. (5.19) describe the interaction correction to the 
total resistance of the system. It is instructive to study 



also the correction to the effective interface resistance de- 
fined in Eq. (4.24). In the limit of strong repulsion, this 
correction is determined by Eq. (5.19a). In the limit of 
weak repulsion, Rt, s — Rr/ sin Q(t) ~ Rr is only a 
small part of Ra at t ^> 1. The A term is also a small 
correction if t <C a/X g . Expanding Eq. ( |5.18| ) in the re- 
gion 1 <C t <C a/Xg and assuming for simplicity that 
A(C) 



A g , one obtains: 



Rt 



err 



Ri 



1 



1 



(5.20) 



where the first two terms come from expansion of 
l/sm9(t). According to Eq. ( |5 20| ), R T , e s(t) first de- 
creases with t up to t ~ (a/A g ) i/4 and then starts to 
increase. This increase is a qualitatively new feature 
appeared due to the interaction. At the upper bound- 



ary of applicability of Eq. ( 5.18 ), at t ~ a/X g , one has 



RT,es(t)/Rr ~ a/Xg ^§> 1. Thus, even within the range 
where the first correction is still small, weak Cooper re- 
pulsion changes the dependence of the effective interface 
resistance which starts to grow with the increase of t. 
The behavior of i?T,eff(i) at t > a/X g will be studied in 



Sec. V C S 



2. Strong repulsion 

If the Cooper interaction is strong in the sense that 

Gt 



X„ > a : 



(5.21) 



then the initial stage of the evolution of u(x) is better 
represented in terms of th e R G equ atio ns for the coeffi- 
cients 71 and 72, cf. Eqs. (3.£) and (5.7): 



^C~ = - A(C)71 ' dc" = -I 71 ' 



(5.22) 



where A(C) is given by Eq. (5.9). Simple calculation gives 



7i(0 
72(C) 



aX 



X g cosh X g ( + Xd sinh A S C 
a 2 /4 



X g cothA„C + Ad' 



(5.23a) 
(5.23b) 



At C 3* Ci = A g , the first harmonics 71 vanishes, while 



the second one saturates at 72 



(i) 



-a 2 /4(A d + A g ). The 



corresponding subgap conductance is then given by 

(5.24) 



Grp 



ATig{X d + X g ) ' 



At the scale C > Ci higher harmonics 7 n >2 are still much 
smaller than 72 once the condition (5.21) is fulfilled, so 
that 
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u(x) — — asin2x, 



2(A d + A 9 ) 



(5.25) 



From now on further evolution of u(x, £) with the in- 
crease of C is given by the solution of the Eulcr equation 
(3.11) with the initial condition ( 5.25| ) [a somewhat sim- 
ilar situation will be discussed in Sec. [VI A |. Change of 
variables u(x, £) = v(y, rf) with y — 2x — it and r\ = 2Q 
makes the solution for v(y, rj) identical (up to the re- 
placement a — » a) to the solution for u{x 1 Q) discussed in 



Sec. Ill B . In particular, the "shock wave" develops in the 
solution for u(x) at the scale C2 = 1/25 = (X g + Xd)/ a 2 
that is much larger than £1 under the condition ( 5.21 ). 

To find the subgap conductance of the system with 
ln(L/d) ^ Ci, one should calculate, according to the gen- 
eral rule ( |4.20|) , the quantity u K (|, C) with ( = 2 \n(L/d). 
This calculation is simplified by the fact that u(%,Q ~ 
at £ 3> Ci- A similar situation had been encountered in 
Sec. [V A where the normal-state conductance was con- 
sidered, cf. Eq. (17). Differentiating Eq. ( j5~8| ) over x at 
x = 7r/2 under the condition w(§) = 0, one gets for the 
Andreev conductance the same Eq. ( |4.7| ), with the initial 
condition ( [5.24 ). As a result, one obtains 



1 

Ga 



1 



G 



(i) 



— 14 

2irg a 



(5.26) 



where the first term is defined by Eq. ( 5.24 ), and the 
second term (which can be of arbitrary magnitude com- 
pared to the first one) is the normal-state resistance of 
the film (in units of h/e 2 ). The first term in Eq. ( 5. 26] ) 
has the meaning of the effective interface resistance de- 
fined in Eq. ( [1.24 ). In the leading order over X g /a, i?T,eff 
is given by l/G^ . Retaining also the next-to-leading 
contribution, one can write Rr,eS as 



n 



A7rg(X g + A d ) C(X d /X g ) 



V9 



O 



371 



(5.27) 



where the function C(Xd/X g ) is of the order of one for 



A rf ~ X g . The second term in the RHS of E q. ( J5.27| ) 
is small compared to both terms in Eq. ( |5.26| ) once the 
conditions (5.21) and ln(L/d) ^> d are fulfilled. 



3. Effect of weak repulsion at large scales 



Here we study the solution of the FRG equation (5. J) 
in the limit of weak repulsion. For the sake of sim- 
plicity we will consider the case of scale-independent 
A(C) = const corresponding to the Finkelstein's fixed 
point A = X g . 

Consider first the qualitative effect of weak repulsion, 
Xg <C a, on the evolution of the function u(x). At small 
t <C a/Xg, it leads to a small reduction of the amplitude 



of the solution compared to the noninteracting case. At 
t ~ a/Xg ^> 1, when this correction becomes of the order 
of the solution itself, u(x) changes its sign and becomes 
negative on some part of the interval < x < it/2. For 
even larger t a/X g , u(f) quickly approaches zero so 
that u(x < |) < and u(x > §) > 0. Then the char- 



acteristics (5.11) with x > it/2 keep going on the right, 
while the characteristics with x < it/ 2 go on the left, 
to the direction of negative x. As a result, shock wave 
will appear also at x = that, however, does not influ- 
ence any physical results given by derivatives of u(x) at 
x = tt/2. The numerical solution for A(£) = A g = a/2 is 
shown in Fig. |t]. 
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FIG. 7. Solutions for Eq. (gj) with A/a = 1/2 for different 
values of t: a) initial stage of evolution; b) solution for large 



To analyze Eq. (gj) in the limit X g <C a, it is conve- 
nient to make the following rescaling: 



s = Xg(, v = U/Xg 



(5.28) 



Then, in terms of the function v(x, s), Eq. (5.8) takes the 
form 

v s + vv x = — 1>(§) sinx, (5.29) 
with the initial condition vq(x) = {a/X g ) smx. 
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FIG. 8. The universal function T(s) = Rt,c«/Rd vs. 
s = (\nL/d)/-Ky/g. 

Since the magnitude of this initial condition is much 
larger than 1, v(x,s) acquires a saw-like behavior at 
s ~ Xg/a -C 1 (i. e., at t ~ 1). Therefore, the details 
of i?o(a;) are irrelevant for the study of the solution at 
large scales, s ~ 1, and one can formally consider the 
initial conditi on vq (x) — Ax with A — > oo. The solu- 
tion for Eq. ( HH ) obtained with such Uo(&) does not 
depend on A g and describes a universal behavior of the 
system for i > 1 in the limit of weak repulsion. The 
corresponding dependence T(s) = Rt, c s{s)/Rd of the 
effective interface resistance obtained by numerical so- 
lution of Eq. ( J5.29 ) and normalized to Rt> is shown in 
Fig. @. The function T(s) has a maximum T = 0.406 at 
s = 1.95, and in the limiting cases it is given by 



T( S ) 



(1-2/tt)*, fors«l; 
1.19/s, fors>l. 



4- Arbitrary A and t 



(5.30) 
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FIG. 9. Dependence of the effective interface resistance 
RT,es(t)/ Rt vs. t — Rd/Rt for diff eren t values of A/a ob- 
tained by numerical solution of Eq. ( |5.8|) for \{Q = const. 
The dashed line corresponds to the noninteracting case, A = 0. 



The most striking feature of Fig. ^ is a nonmonotonous 
dependence of i?T,cff(i), especially pronounced for weak 
repulsion, X g -C a. In this limit, the effective interface 
resistance significantly exceeds its noninterating value, 
Rt, at large scales, t > (et/Ag) 1 / 2 . Such a nonmono- 
tonuos dependence arises even within applicability of the 
first- order correction ( 5.1 8|) as discussed in the end of 
Sec. V C 1 . Another important feature of Fig. || is that 
the limits A — ► and Ro — > oo do not commute. Indeed, 
for any small but finite A g , Rt, c s (t)/Rr will eventually 
(though, at very large t) deviate from the noninteracting 
dependence (dashed line in Fig. ||) and become large. 



We will sec in Sec. VI that such a nonmonotonous 



dependence of RT, e s(t) will manifest itself in a non- 
monotonous temperature and voltage behavior of the 
subgap conductance. 



Solution of the FRG equation (5.S) for arbitrary A and 
t should be obtained numerically. The effective interface 
resistance Rt.cS normalized to the tunneling resistance 
Rt as a fu nction of t = Rd /Rt is plotted in Fig. || (as in 



Sec. V C 3, we consider A(£) = \ g = const). The dashed 



line shows Rt.cS / Rt = 1/ sin 0(f) for the noninteracting 



case [cf. Eq. fl4.23|)]. 

For the case of strong repulsion, \ g 3> a, i?T,eff(f) 
very quickly (at t ~ a/X g <C 1) reaches its asymptotic 
value and saturates at i?r,cff(f = oo) ~ (2A s /a)i?r, cf. 
Eq. (5.27). The limiting value Rt, c s{oo) decreases with 
the decrease of X g /a up to X g /a ~ 1. For small X g /a, 
corresponding to the case of weak repulsion, it starts to 
grow ag ain w ith i?<r,eff(oo)/i?T ~ 1.19a/A s , according 
to Eq. (5.30). Note that in this limit i?T,eff(i) reaches 



its asymptotic value at large scale t ~ &/X g . The de- 
pendence of i?T,cff (oo) as a function of X g /a is shown in 
Fig.0. 




FIG. 10. The limiting value Rr,es{t = oo) /Rt vs. A/a (for 
X(C) = const). 
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FIG. 11. The effective interface resistance Rr tC g normal- 
ized to its limiting value at G^, 1 = oo (see Fig. hd) as a func- 
tion of y/g/Go = A 9 £/2 for various values of Gr/y/g- The 
plot corresponds to the Finkelstein's fixed point A = A 9 = 

In this section we assumed that A(£) had already 
reached its fixed-point value X g = l/ftiry/g), so that 
\ g /a = Thus, the limits of strong and weak 

interaction are separated by Gt ~ \f§- For Gt <C ■Jg, 
interaction is strong, and for Gt 3> \fg (and, in par- 
ticular, for the case of a completely transparent inter- 
face, Gt = oo), interaction is weak. It is then natu- 
ral to measure all dimensionless conductances in units of 
^fg. The effective interface resistance Rt, c s normalized 
to its limiting value i?T,eff(oo) at large spatial scales (cf. 
Fig. |l^) as a function of y/g/Go = s/2 = X g (/2 is shown 
in Fig. [ll], with different curves corresponding to dif- 
ferent values of Gt/ ' \fg~- The asymptotic curve (dashed 
line) for the transparent interface, Gt — > oo, is given by 
R TcS = T(2^/g(e 2 R D /h))R D , cf. Eq. ( pHSC) ). In this 
limit, for Rtj ^> h/e 2 y/g one has 



R - °' 6 

ft Vg 



— *4 

2ng a 



(5.31) 



D. Noise 

In this section we will analyze the effect of interac- 
tion on the noise of N-S current. As shown in Sec. IV C , 
the current-current correlator is determined by the noise 
function Ps{t), see Eqs. ( 4.25 ) and ( 4.26 ). Qualitatively, 
Ps can be estimated by comparing the effective tunnel- 
ing resistance, Rt^s, with the diffusive resistance, Rd- 
If i?r,eff S> Rd then Ps « 3, while in the opposite case 
(Rt,c« < Rd) Ps ~ I- As in Sec. |VC3| , we will as- 
sume here the A(£) had already reached its fixed point 
A = X g . Then Ps becomes a function of two parame- 
ters: t = Rjj/Rt and X g /a — 2^/g/GT- Summarizing 
results of Sec. |V C| , one can sketch the boundaries be- 
tween regions with Ps = 3 and Ps — 1 on the plane 
(logt, log(G T /V5)), see Fig. e 




FIG. 12. Schematic map of the noise coefficient Ps as a 
function of the tunnel and diffusive resistances in the pres- 
ence of Cooper repulsion A = A s . Dark area refers to the 
tunnel limit Ps fa 3, whereas light regions correspond to the 
diffusive regime with Ps ~ 1. 



In the case of strong repulsion, Gt <§C \/g, the func- 
tion u(x) very quickly ( at t ~ Gr/y/g) reduces to the 
second harmonic, Eq. ( 5.25j ), and Rt, c s saturates at 



i?T,cff = R { a ] = (4J g/G T )R T [cf. Eq. J5.24D 1. As ex- 
plained in Sec. V C 2[ further evolution of u(x) is analo- 
gous (after a proper rescaling and shift of variables) to the 
evoluti on of u(x) for the noninteracting case considered 



in Sec. Ill E| . The property u(§) = makes the calcula- 



tion of the noise function Ps analogous to the calculation 
of the function P/v for the case of the normal island, see 
Sec. fVA[ So, one obtains P s nt = P n (R d /Ra ) )- Thus, 
in the limit Gt <SC \fg, the crossover between the tunnel 
(Ps = 3) and diffusive (Ps = 1) character of noise is 
shifted to t ~ y/g/G T > 1, see Fig. ||. 

In the limit of weak repulsion, Gt 3> y/g, the situation 
is more interesting. For t ~ 1, interaction corrections can 
be neglec ted a nd Ps is given by the noninteracting ex- 
pression ( 4.2S ). So, at t ~ 1, Ps decreases from 3 to 
1, the corresponding boundary being shown in Fig. 
Later, at t ~ Gt/ ' Jg 3> 1 (when Rd ~ h/e 2 y/g) inter- 
action corrections become relevant. In this region, i?T,eff 
is of the order of Rd (cf. Fig. |||), and one may antici- 
pate that Ps will deviate from 1. For even larger t when 
resistance is dominated by the diffusive conductor, Ps 
will eventually reduce down to 1. This crossover region 
is marked in Fig. 12 by the dashed lines. 



The behavior of Pg in the crossover region t ~ 
Gt / ' \fg 3> 1 can be obtained by numerical solution of 



Eq. (5.29). The resulting Ps(s) is plotted as a function of 
s = X g C = (hxL/^/ny/g = 2(e 2 /h)R D yfg [cf. Eq. " 
in Fig. Ps{s) has a mimimum Ps = 0.99 at s 
and a maximum Ps = 1.28 at s = 3.25. 



5.28)] 



0.40 
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FIG. 13. Noise function Ps(s) vs. s = for the case 

of weak interaction, Gt 3> ^f§- 

The negative derivative of Ps(s) for s <C 1 canbe ob- 
tained analytically. Seeking the solution of Eq. ( 5.2S ) as 
a series over s, one finds 

x 7r(l — cos x) „ . „ 

= ^— L + 0(s), (5.32) 

s 2a; 

where the first term is the usual saw function, while the 
other terms are due t o the RHS of Eq. ( |5.29 ). Now, with 
the help of Eq. ( 4.26 ) one gets for s«l 



Ps(s) = 1 



12 

^2 



24 \ 

-~ s = 1 -0.058 s. (5.33) 

7T J I 



VI. TEMPERATURE, VOLTAGE AND 
MAGNETIC FIELD EFFECTS 

A. General analysis 

In this section we consider the case of a superconduc- 
tive island in the situation when temperature and/or 
voltage are high compared to the Thouless energy, or 
the perpendicular magnetic length Ijj = \J<&q/H is 
shorter than the system size L. Then Cooperon coher- 
ence is destroyed at the energy scale 51* > -Erh, where 
Q* = max(T, eV, eDH/c). At the same time we will as- 
sume that frequency u> of external voltage (or current) is 
smaller than the Thouless energy. 

According to the general approach explained in 
Sec. Ill C , we have to determine an effective Proxim- 



ity Action via the RG procedure accomplished down 
to the energy scale E^h- Now the whole energy inter- 
val of the RG, ojd > > -ETh, can be divided into 
two regions with different FRG equations. In the re- 
gion u>d > Q > f2*, rcnormalization of the action (3.1) is 



described by Eq. (3.11) derived in Sec. [II for the non- 
interacting case, or its analogue (5.8) in the presence of 
interaction. In the region f2* > Q > i/rh, Cooperons are 
suppressed while diffuson modes still contribute to the 
renormalization of the action (3.1). Thus, at large scales, 



il < fi*, the charges 7„ depend on two RG "time" argu- 
ments, C and C* = ln(wd/ri*). For £ > (*, the bulk matr ix 



Q is purely diagonal in the Nambu space [cf. Eq. (2.3)]. 
On the other hand, according to Eq. ( |2.8| ), the island's 
matrix Qs is off-diagonal in the Nambu space, due to the 
Andreev nature of the subgap tunneling across the inter- 
face. As a result, tr(QsQ) n with odd n vanishes, and so 
7„ 's with odd n do not contribute to the multicharge RG 



equation (3.8) at scales C > C*- 

To describe evolution of 7„'s at £ > one should 
modify the derivation presented in Sec. Ill B for the case 
of even n. The difference is that now only the first line 



of Eq. (A2) corresponding to diffuson pairing contributes 
to the average ASk, n = i^int.fcSint,™)- But for even k 
and n, a nd at £ > C* the matrices A and B defined in 
Eq. (3^) commute with r z . Then the zero-energy limit 
of the first line of Eq. (A2) reproduces Eq. ( |A3| ) used 
in the derivation of Eq. ([3.6]). Hence, we conclude that 
evolution of {72m} at C > C* is described by the same 
Eq. (3^), with all j2 m +i se t to zero. 

To describe evolution of the set {j n } at £ > £*, we find 
it convenient to introduce Fourier transformation with 
respect to even n: 



u(x) 



1, 



00 

2n-f2n sin 2nx — — [u(x) — u(ir — x)]. (6.1) 

n=l 



Then the FRG equation for the function u{x) acquires 
the form of the Euler equation: 



0. 



(6.2) 



The function u(x, C*) that determines physical quanti- 
ties (see below) is then given by the solution of Eq. (6.2) 
with the initial condition 



1 



u{xX*;C*) = oKa;,C*) - u ( n - x,C*)], 



(6.3) 



where u(x,(*) is the solution for the FRG equation (5.8) 
at C* = m (W(i/^*)- We want to emphasize that the re- 
duction u(x) — > u(x) at C = C* describes the crossover 
from the f2* <C -Erh to £1* ^> i?Th regimes only with log- 
arithmic accuracy. The number in the correct cutoff of 
logarithm is beyond the RG precision. 

Below we apply the described scheme to the calcula- 
tion of physical quantities. To determine the Andreev 
conductance at large 11* > £/rh in terms of the function 



u{x) we should use a generalized version of Eq. (4.18). 
Namely, we should take into account that for e > i?Th, 
parameters 72« depend on the energy argument e running 
under the trace [cf. Eq. (B3)], since it is just the value of 
e which determines the Cooperon coherence scale. Thus, 
72n(e) should be put under the sign of Tr in Eq. (4.18). 



The trace operator contains an integral over e those main 
contribution comes from e *~ fi*. As a result, we obtain 
the following expression for the (nonlinear) current: 
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I(V) = —G A (Sl,)V, 



(6.4) 



B. Noninteracting case 



where the value of G A (fl*) is determined to logarithmic 
accuracy as 



Ga(O0 = G A (t,U) = 4ngu x (%, C;C*)- 



(6.5) 



Calculation of u x (f) is ve ry simple since u(%) = 0. 
Hence, similarly to Eq. (4/7) one gets 



d(l/G A ) 



1 



(6.6) 



Eqs. (6.5), ( p.q ) can be naturally interpreted with the 
help of the effective interface resistance Rt.cS introduced 
in Eq. (4.24). An important property of this quantity 
is that it is formed by Cooperons only which are taken 
into account by the FRG equation for u(x). At scales 
t > t* = a£*, Cooperon coherence is lost, and i?T,eff(i) 
saturates at the constant level i?T,off(i*) becoming inde- 
pendent of Rd any more. As a result, 



R A (t, t*) — Rd + RT,eff(t*), 



(6.7) 



where Rd — Rrt = \\T.{L/d)/2ira is the energy- 
independent resistance of the normal film. 

Consideration that lead to Eq. (|6.5|) can be generalized 
to higher correlators of current as well. To logarithmic 
accuracy, the noise power is given by the zero-energy 
expressions ( 4.25| ) and ( 4.26 ) with u(x, £) replaced by 
u(x, C;C*)- The quantities U x (%) and U XX x(%) en terin g 
this expression can be calculated similarly to Eq. (4.14): 



^x { n J C) C* ) 



u»(f,C) 



L xxx \ 2 



( 2 ' Cs C*) 



i + (C-C*K(f,C*)' 

Uxxx ( 2 ) C* ) 



(6.8) 
(6.9) 



[i + (C-CK(i>C)] 4 ' 

As a result, the current-current correlation function can 



be written in the form similar to Eq. (4.25) 

e 2 G A (t,U) 



:m 



{(3 -P s (t, *.))*(<") 



~Ps(t,U) - 2eV) + 9(w + 2eV)]}, (6.10) 



where G A (t, t*) is given by Eq. (6.7), and the noise func- 
tion Pg(t,t*) depends now on two RG "time" variables. 
Using Eqs. (gjoh, flOg 1 ), O) and 0, one obtains 



G^(t, t*) 
G^(t.,t.) 



(6.11) 



Here both GA(£*,t*) - G A (t*) and Pg(t*) are given by 
the zero-energy results at t = i*. 



_Z. ylnrfreeii conductance 



To calculate the Andreev conductance we substitute 



R TtCS (U) = R T /smQ(Q into Eq. (4.24) 



l + tsin©(£» 



l + tsin0(£» 



(6.12) 



that can be obtained from the zero-energy result (4.23) 
by replacement 0(f) — > 0(i*). 

Below we apply Eq. ( 6.12| ) to the analysis of two spe- 
cific examples. 

We start from the case of the linear conductance as 
a function of temperature. The corresponding curves 
G A (T) for several values of the ratio t — Gt/Gd are 
presented in Fig. [li]. In the limit Gt ~> Gd any depen- 
dence on T disappears and N-S conductance is equal to 
the diffusive conductance Gd- This result is in disagree- 
ment with calculations |33],[l7) which predict, for an ideal 
N-S structure in the ID geometry, a conductance maxi- 
mum at T ~ ETh with the relative magnitude about 10% 
(the so-called finite-bias anomaly). This discrepancy is 
due to the limited precision of our calculation scheme, 
which is equivalent to the summation of the main loga- 
rithmic terms. Finite-bias anomaly is due to the energy- 
dependence of the effective diffusion constant D(e) at 
e ~ i?Th, and this effect is beyond the main logarithmic 
approximation. As one can see from Fig. [l4], the zero- 
bias anomaly is stronger than the finite-bias anomaly for 
small enough values of t. 




FIG. 14. The linear Andreev conductance G A (T) normal- 
ized to its zero-temperature value vs. T/Etii for different val- 
ues of the ratio t = Gt/Gd- The plot corresponds to the 
geometry shown in Fig. [j] with L/d — 20, i. e., Ud/E^ — 400. 
In the intermediate region T ~ _&rh (sketched by dashed line) , 
corrections are nonlogarithmic and cannot be taken into ac- 
count within the RG approach. 
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Consider now the nonlinear sub-gap conductance at 
high temperature T ^> Et\i and arbitrary relation be- 
tween T and eV. At eV «Twe are back to the linear 
conductance case with = = {Gx/^g) ba(wd/T), 
whereas at large voltage e7 > T we have i* = ty = 
(G T /^g)\n(u d /eV). 

To find the behavior of Ga(T, V) in the crossover re- 
gion eV ~ T » i?Th we need to improve the logarith- 
mic accuracy of Eq. ( |6.12| ). To this end, we perform a 
more accurate calculation of the energy integral under 
the trace in the generalized version of Eq. ( 4.1 8j ) taking 
into account the energy dependence of j n (E): 



Ga(T,V) 
Gd 



1 



dE 



tanh ■ 



2eVj Q 

tsm<d(t E ) 
X l + tsmQ{t E )' 



E+ 
2T 



tanh ■ 



E_ 
2T 



(6.13) 



where E± = E±eV, and t E = (G T /4ng) hx(u d /E). The 
logarithmic factor \a.{uJd/E) comes from the integration 
over 2D Cooperon modes, J 1 d 2 q/ (Dq 2 ± 2iE), which 
determines the Cooperon amplitude at coinciding point 
(i. e., the probability of return). In the presence of trans- 
verse magnetic field, Cooperon modes are quantized, so 
integration over momenta is substituted by the summa- 
tion over Landau levels, see, e. g., Ref J34|. Then the 
effect of magnetic field upon the subgap conductance can 
be accounted by the replacement of t E in Eq. (3.12) by 




2eV/T 

FIG. 15. The nonlinear Andreev conductance Ga{T,V) 
normalized to the linear conductance Ga{T, 0) vs. the ratio 
2eV/T for T/E Th = 10. As in Fig. [L4| u d /E T h = 400. Dif- 
ferent curves correspond to different ratio t = Gt/Gd'- 0.1 
(bottom), 0.5, 1, 2, 10 (top). 

For arbitrary t, the behavior of G a in the crossover 
region should be obtained by numerical integration in 
Eq. (6.13). The corresponding plots of Ga(T,V) as a 
function of 2eV/T for different values of t are presented 
in Fig. |l^. Note, finally, that the difference between the 
nonlinear conductance Ga(V) = Ia/V and the differen- 
tial conductance dlA/dV should be neglected in the main 
logarithmic approximation. 



2. Current fluctuations 



t E (H) 



Gt 
Airg 



In 



he 



2eHd 2 



Retp 



. Ec \ 
' 2eDH J 



(6.14) 



where ip( x ) — d\nT(x)/dx is the digamma function. 

In the limiting case of weak interface transparency, 
t <C 1, the Andreev conductance is given by 



G A (T,V,H) 



eV ' 
2ujd 
irT 



b 1, for eV > (T, ^); 
+ 7, for T » (eV, ^); 



In ^+ 7 , for^»(T,eV); 



eDH 



(6.15) 



where 7 = 0.577 ... is the Euler constant. Comparing 
these asymptotics, we conclude that crossover from the 
voltage to temperature-dominated effective interface re- 
sistance occurs at 2el^ « 7re 1_7 T ss 4.8 T. Similarly, 
crossover from the temperature- to magnetic-field domi- 
nated resistance occurs at H « TrcT/eD. 



In order to calculate the noise function Ps(t,t*) at 
fi* > E T h we substitute Eqs. ( 4.2S ) and ( 6.12j ) into 
Eq. ( J6 . 1 1[ ) . As a result, we obtain 



1 + 6* tan 6* +36* cote* 
Ps(M * ) = 1+ 2(1 + 6* tan6*)(l+isin6*)3< ((U6) 

where 6 = Q(t), 6* = 6(i*), and G^fi*) is given by 
Eq. (ph. 
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FIG. 16. Noise function Ps(t,tt) vs. S7,/i?Th for different 
values of t — Gt/Gd- As in Fig. u4, the ratio Lo d /E T h = 400. 
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Eqs. (6.10) and ( |6.16| ) determine fluctuations of N-S 
current as a function of temperature, voltage, magnetic 
field and frequency (uj <C £rh) at arbitrary values of 
the ratio t — Gt/Gd- Dependence of noise upon T, 
V and ui comes in two different ways: via the func- 
tions *(w ± 2eV) = Slcoth((w ± 2eV)/2T) and via 
the Cooperon cutoff scale S7* = max(T, eV, u>, eDH/c) 
which determines t* = (Gt /47rg) ln^^/Sl*). In the limit 
t — > oo, the function Ps(i,i*)/3 approaches 1/3 and the 
expression for noise reduces to the standard form [jl9| for 
purely diffusive N-S junction at SI* < E^h ]§H|- How- 
ever, the range of t where such a universal behavior sets 
in, does depend upon SI*: the increase of SI* leads to 
the decrease of the proximity angle 0*, which, in turn, 
increases the factor Ps (i , i* ) . A number of curves charac- 
terizing the behavior of Ps{t, <*) as a function of Sl*/£>rh 
at different values of t is presented in Fig. [L6L 



C. Effects of interaction 

This section contains the main application of our the- 
ory since variation of temperature, voltage, or magnetic 
field is a natural tool to study system properties. Here all 
effects discussed in the body of the paper come into play 
altogether. Therefore, we will repeat briefly the main 
concepts the reader could gain from the above discus- 
sion. 

At scales smaller than the Cooperon coherence length 
yjD/£l*, i. e., at C < C* = ln(wd/Sl*), the system 
is described by the function u( x, C ) which evolves ac- 
cording to the FRG equation ( |5.8| ). At larger scales, 
C > C*j onc should introduce a "two-time" function 
u(x, C ; C* )- It is obtained from u(x, Q by the reduc- 
tion ( |6.3|) and obeys the FRG equation ( |6.2|) . The con- 
ductance and noise power are given by Eqs. ( 4.20| ) and 
(4.26), with u(ir,£) being substituted by u(x, C;C*)i and 
( = ln(ujd/ Et\i) = 2ln(L/d). The fl* dependence of the 
conductance can be easily expressed with the help of the 
effective interface resistance according to Eq. (6.7). 

In the noninteracting approximation the effect of large 
S7* i?Th was to decrease the Andreev conductance 
compared to the zero-energy limit, SI* <C P-rii, due to the 
increase of i?T,cff(i*)- The drastic change introduced by 
the Cooper interaction is that i?r,eff(i*) shown in Fig. ^ 
is no longer a monotonous function of t*. Therefore, de- 
pending on the relation between parameters of the prob- 
lem, the sub-gap conductance may either decrease or in- 
crease with the increase of SI*. This unusual enhance- 
ment of conductivity with the increase of the decoher- 
ence energy scale SI* is most pronounced in the limit 
of weak repulsion, Gt ^> 2^/g (we again assume that 

In 



VC1 



A(C) had reached the Finkclstcin's fixe d poin t X g ) 
this case, according to the result of Sec. 
decreases with t* at t* <c {Gr/^/g) 1 ^ 
with U at U > (G T /V<?) 1/4 , cf. Eq. 



total conductance decreases with the growth of SI* for 
SI* 3> Sl cr and increases for SI* <C Sl cr , where Sl cr — 
cu d exp(-cg 7 ^/G T /4 ), and c = 27r 7 / 4 (7r - 2)~ 1 / 4 « 14.3. 
Since th e to tal resistance is the sum of Rt, c s and Rrj, 
cf. Eq. (|3), the magnitude of the effect is determined 
by the ratio Rt.cs/Rd which, according to Fig. ||, has 
a maximum at t ~ Gt / ^fg- An example of such a non- 
monotonous dependence of Ga(Q*) is shown in Fig. [l?]. 
The curves differ by the ratio t = Rjy / Rt and correspond 
to G T = 10^5 (i- e., \ g /a = 0.2). 




FIG. 17. Dependence of the Andreev conductance Ga(Q*) 
(normalized to the zero-energy value) on the ratio 
for different values of t. For all plots, o^/Stii = 400, A = X g , 
and Gt = W^/g- 

In the opposite case of strong repulsion, Gt <C 2^fg, 
the SI* dependence of Ga is absent for i* ^> GT/y/g 
when RT, e s(t*) is i* independent. The reason is that 
strong repulsion makes Cooperons ineffective at scales 
C ^> Ci = 1/A g when all 7„'s with odd n vanish, see 

u(x) at the 



Sec. |VC 2[ Therefore, the reduction u(x) 
scale £* leaves the function u(x) intact indicating that 
physical results are SI* independent. They become S7* 
dependent at relatively large scales when £* < d, i. e., 
at ln(w<j/Sl*) < 2iT^/g. The principal effect of SI* is then 
to decrease Ga with the increase of SI* . 



RtMQ 

and increases 
( pTlil ). Hence, the 




FIG. 18. Noise function Ps(s,s*) vs. s, for different values 
of s. The dashed line shows Ps(s*, s„) = Pg(s*). 
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Finally, we will dwell on the fi* dependence of the 
noise power. According to Eq. ( |6.10 ), the latter is de- 
scribed by the function -Ps (£;£*) which now depends on 
two RG "time" arguments. Eq. ( |6.11 ) relates Pg(t, £*) to 
the zero-energy expressions for GU(i*) and Pg(t*). 

We start from the case of the weak interaction. The 
crossover from the tunnel (Pg = 3) to diffusive (Pg = 1) 
character of noise at t, £* ~ 1 is well described within 
the noninteracting approximation investigated above. To 
study the noise function in the vicinity of the interaction- 
induced peak in Pg at t, t* ~ Gt/ y/g S> 1 (cf. Fig. |l^) we 
substitute the functions T(s) a nd Pg (s) shown in Figs. ^ 
and [l3|, respectively, into Eq. fl6.ll ). The resulting de- 
pendence of Pg(s,s*) as a function of s* is plotted in 
Fig. [l8]for different values of s. We remind that Ps(s, s») 
is defined at s* < s and reduces to the zero-energy result 
at coincident "times": Ps(s*,s*) = Pg(s*). The latter 
function is sketched in Fig. [l8| by the dashed line. Tak- 
ing into account that s* = A 9 £* = (l/2iTy/g) Infad/fil*), 
one obtains from Fig. [l^ that Pg decreases with the in- 
crease of the Cooperon decoherence energy scale f2* (at 
t, t„ ^> 1), as if the system becomes more diffusive. This 
trend is opposite to what one has in the noninteracting 
case when the increase of fi* drives the system toward 
the tunnel limit, thus, increasing Pg. Note also that, 
contrary to Ps(s) (cf. Fig. 13), Pg(s, s*) shown in Fig. [l8| 
is a monotonous function of at a fixed s. 

In the limit of strong repulsion, Gt <C y/g, the zero- 
energy noise function Pg(t) exhibits a crossover from the 
tunnel to diffusive regimes at t ~ y/g/Gr ^S> 1. Neverthe- 
less, Pg(t,t*) remains t* independent upto much smaller 
i* ~ Gx/y/g <C 1 corresponding to relatively large en- 
ergy scales, \n(u)d/£l*) — ^yfg- This effect is of the 
same origin as the above-mentioned Q*. independence of 
Ga in the limit of strong repulsion. 



VII. CONCLUSIONS 



We have shown in this paper that electron transport 
through mesoscopic N-I-S structures can be fully de- 
scribed in terms of the Keldysh-space Proximity Action 
S P iox[Qs,Qn], as defined in Eq. fl3.14|) . This action is a 
functional of two matrices, Qg and Qn, corresponding to 
the superconductive and normal leads. Throughout the 
paper, we choose the gauge with the normal lead being 
in equilibrium, with Qn — A. The superconductive ter- 
minal is characterized by the matrix Qg which contains 
both the classical (<fi) and quantum (<f2) components of 
the order parameter phase. Physical response and cor- 
relation functions of any order can be determined from 
the Proximity Action by calculating the derivatives with 
respect to the quantum component at a given value of 
the classical component of the phase, as explained in the 



functional contains information about full charge transfer 
statistics (cf. Refs. @||). 

The Proximity Action functional is known once the 
set of "charges" "f n or, equivalently, the periodic func- 
tion u(x), Eq. (3.1C), is specified. It was explained 



in Sec. Ill C that the function u(x) is directly related 
to the generating functional for the transmission coeffi- 
cients, introduced in Ref. B. Actually, the Proximity 
Action approach bears an obvious analogy with the scat- 
tering matrix approach as both describe transport 
properties in terms of the characteristics of the termi- 
nals (stationary-state Green functions of the terminals 
Qg.N in the former versus asymptotic scattering states 
in the latter approach). The Proximity Action method 
is well-suited for the treatment of interaction effects in 
the contact region, the task which is out-of-reach for the 
standard scattering-matrix technique. 

Actual calculation of the function u(x) that determines 
the Proximity Action is accomplished by the Functional 
Renormalization Group (FRG) method. In the sim- 
plest case (no interaction in the N region, and all rel- 
evant energies are much below the Thouless energy scale 
-Erh = fiD/L 2 ) the FRG equation reduces to the Eu- 



ler equation (3.11) for ID motion of a compressible gas. 
Although we derived this equation for the normal con- 



beginning of Sec. [V. In principle, the Proximity Action 



ductor of 2D geometry, its solution, Eq. fl3.13| ), is appli- 
cable (being expressed in terms of the ratio t = Gt/Gd 
of the tunneling to diffusive conductances) to any co- 
herent conductor. If higher energy scales £1 > i?xh or 
electron-electron interactions are involved the applicabil- 
ity of our FRG procedure is limited to 2D diffusive con- 
ductors. The generalized FRG equation that takes into 
account interaction constant A in the Cooper channel of 
the N conductor is given by Eq. fl5.S|). We derive this 
equation and analyze its solution in Sec. analytically 
in two limiting cases of weak interaction of arbitrary sign 
and of strong repulsion. We also provide the results of nu- 
merical solution of Eq. (5.£) in the intermediate region of 
moderate repulsion. The main feature of the interaction 
effect in 2D is that it leads to the infrared-growing cor- 
rections of the order of A \n(L/d), that should be summed 
up non-perturbatively for the system of sufficiently large 
size L, for any value of the ratio Gt/Gd- This is in con- 
trast with the ID case (as treated in Ref. [jl7j for the case 
of the fully transparent interface, G^ 1 = 0), where the 
interaction correction to the resistance was found to be 
small, of the order of A itself. 

Physical consequences of Cooper repulsion in a 2D 
proximity system depend on relation between A and the 
ratio Gr/^g- For clarity, we assumed that interac- 
tion constant had reached the Finkelstein's fixed point 
A = l/27Ty / g. Thus, strong repulsion corresponds to a rel- 
atively weak tunneling conductance Gt <C 2y/g, and vice 
versa. It is convenient to characterize N-I-S transport in 
terms of an effective interface resistance Rt^b, so that 
the total resistance is equal to Rd + Rt, c b- If Gt <C 1\fg~, 
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then i?T,eff ~ (4y / 5 '/ 'Gt)Rt 3> Rt- This asymptotic re- 
sult is valid at t > Gt/ yfg- At longer length scales, such 
that ln(L/d) > n-y/g, Cooper-channel correlations are ir- 
relevant and the form of the Proximity Action coincides 
with that of an N-I-N structure, up to the replacements 
of charge quantum and the tunneling resistance: e — > 2e 
and Rt — ► Rt.cS- As a particular example of this prop- 



erty, in Sec. VD we calculated the N-I-S noise intensity. 
In the opposite case, Gt ^> 2y/g, the evolution of Rt.cB 



with the growth of t is rather unexpected (cf. Fig. 11): 
first it decreases with t, reaches a minimum value of the 
order of Rt at t ~ (Gt/ y/g) 1 ^ 4 , and then grows up to 
the asymptotic value 0.6 H/e 2 ^/g at t > Gr/y/g- Non- 
monotonous behavior with a maximum of relative height 
about 30% is demonstrated by the noise power as well, 



cf. Sec. |VD 

At finite temperature, voltage or magnetic field 
Cooperon coherence is destroyed at the time of the order 
of 7i/f2*, where 17* = max(T, eV, eDH/c). We calculated 
the Andreev conductance and shot noise power in the 
presence of these phase-breaking effects in Sec. VI B for 
a non-interacting normal conductor. The role of phase- 
breaking effects in the presence of repulsion in the normal 
conductor is discussed in Sec. VI Cj . In the case of highly 
transparent interface, Gt ^> 2yfg, and at Rd 3> Rt, an 
unusual "finite-bias" maximum of the Andreev conduc- 
tance Ga(Q*) is predicted, cf. Fig. [l7| which is of the con- 
sequence of non- monotonous behavior of Rt off (t) men- 
tioned above. Contrary to the well-known |33|,[l7| finite- 
bias anomaly at (T, eV) ~ i?Th in the non-interacting 
case, this new anomaly is expected at much larger value 
of f2*/-Erh- Moreover, it should be seen as a function 
of magnetic field as well. In simple terms the origin of 
this new effect can be understood as follows: repulsion in 
the normal metal produces a superconductive "gap func- 
tion" in the normal conductor, Ajy, with the negative 
(compared to As in a superconductor) sign. Due to its 
opposite sign, A^r decreases the conductance of the struc- 
ture, therefore any decoherence that reduces A at leads to 
the increase of the conductance. 

We did not consider in this paper the weak localization 



(WL) effect and the Coulomb zero-bias anomaly (ZBA). 
The reason to neglect them, as explained in the Intro- 
duction, is that relative magnitudes of these effects scale 
as <7 — 1 1 In 1 . However, both of these effects have high- 
frequency cutoff at the inverse elastic scattering time 1 /r, 
and their contributions to the tunneling and diffusive 
conductances from the energy scales cud < ^ < 1/t are 
not necessarily small. Since logarithmic corrections to 
the Proximity Action we have studied are coming from 
much lower energy scales ft < uod, these effects can be 
considered separately: high-frequency fluctuations can 
be accounted via the replacement of the bare conduc- 
tances Gt and g by their values corrected by the WL 
and ZBA effects at the energy scales uod < fi < 1/t, 
provided that the renormalized value of g is still large. 

In this paper we consider the low-frequency case lu <C 
i?Th only. This limitation is due to technical complica- 
tions: at higher frequencies the effects of dynamic screen- 
ing should be taken into account, which makes the ex- 
pression for the Proximity Action much more involved. 
However, it is possible to show that the results obtained 
for the Andreev conductance are still valid for frequen- 
cies up to much higher frequency scale w max [defined af- 
ter Eq. fl4.9p ] at which the electroneutrality of the system 
sets in. Experimentally, N-S noise in the presence of high- 
frequency radiation (to ^> £rh) was studied in Ref. fl36| . 
We leave theoretical consideration of this problem for the 
future. 
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APPENDIX A: CONTRACTION RULES 

The contraction rule for averaging over W is given by [Q 

where 

tr [A 12 A 2 B 2 iA 1 + A 12 (A ) 2 B 21 (A ) 1 - A 12 t z B 21 t z - A 12 B 21 ] [Dq 2 + - e 2 )(A )i] 



X 



(Dq 2 ) 2 + (ei - e 2 ) 2 

tr [A 12 A 2 B 21 A! - Ai 2 (Ao)2 J B 2 i(A )i + A 12 t z B 21 t z - A 12 B 21 ] [Dq 2 + i(e x + e 2 )(A )x] 



(Dq 2 ) 2 + ( Cl + e 2 ; 



2 



(Al) 



(A2) 
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Here A 12 = A(q, e 1 ,e 2 ), B 21 = B(-q,e 2 ,e 1 ), A m = A(e m ), and (A ) m = A (e m ). The first (second) line of Eq. ([A|) 
corresponds to diffuson (Cooperon) pairing. 

In the zero-energy limit, Dq 2 3> ei, e 2 , this expression can be simplified as 



(TtAW ■ Tr BW) = ± f tr ^ A ^ A * ~ A M 



■Kg 



(A3) 



APPENDIX B: CONTRIBUTION OF A SINGLE 
TERM FROM THE MULTICHARGE ACTION TO 
THE CONDUCTANCE AND NOISE 

Here we calculate the quantity which describes the con- 
tribution of a single term from the action (3.16) to the 
DC conductance of the system in the zero-energy limit 
(A, < £ T h): 



— Tr^( e ^A e"^Ao) r 



= - n 2 eV. (B4) 

¥2=0 7T 



In a similar way one can obtain 



— TrK(e^ /2 A e-^/ 2 A )" 
dtp 2 



= - n 2 eV. (B5) 



— Tr^(e^A e-^A y 
dtp 2 



V2=0 



(Bl) 



Now we turn to another expression emerging in calcu- 
lation of the current-current correlator from the action 



with ipi obeying the Josephson relation ipi = 2eV. The 
derivative can act either on e tv or on e~ %v so that (Bl) 
reduces to 

f°° de 

in / — tr K a x [M n (V) - M n {-V% (B2) 
where M is diagonal in the energy space matrix: 
M ee (V) = A (e+)A (e_) = 



L n {t,t';V) 



SMt) 6<p 2 (t') 



ip 2 =0 

(B6) 



1 2[F(e-)-F(e+)] 
1 



and e± — e±eV. Calculating M°, integrating over e and 
tracing with a x one finds that both terms in Eq. (B2) 
yield the same contributions, and 



Here there are several possibilities depending on where 
the derivatives act on. First of all we note that they can- 
not act on the same exponent e ltp (or e~" ltp ). Indeed, 
in this case Eq. (jB6j) would reduce to TrM n (V) cx Tr 1 
, that would give zero according to the rules of the Keldysh 
k a- model [^T|. Therefore there are four different contri- 
(B3) butions (coinciding time indices t, i! are omitted): 

L n (V) = K(V) + L a n (-V) + L b n (V) + L b n (-V), (B7) 

where 



J 



L»(t, £■ V) = n ]T Tt k e 2leVt A M p (—V) a x e- 2leVt ' A M q (V) a x 



p,q>o 

p+q=n—l 



L b n (t,t';V) = -n J2 Tr K e 2ieVt A e- i ^A M r (V)a x e 



2ieVt' 



A e~^ A M s (V)a x 



(B8) 



r,s>0 
r+s=n-2 



Performing Fourier transformation to the frequency domain one obtains 

/OO J 
— tV K (A (-10 ^)e+u>+2eV (A M*(V) O x ) €l 
-oo M 



p,q>0 



/°° de 
— tr K A (e + uj + 2eV) (A M r (V) a x ) e+w A (e + 2eV) A M s (V) a x ) e . 
-CO 



(B9) 



s>0 
r+»=n-2 



Evaluating matrix exponents and computing traces we obtain 

/OO 
[(p + l)F(e + lo + 2eV) - pF{e + uj)} [(q + l)F(e) - qF{e + 2eV)] 
CO 



p,q>o 

p + q = n, — l 



/CO 
[F(e + u) - F(e + u + 2eV)] [F(e) - F(e + 2eV)] 
-00 



s>0 
r+s=n-2 



27' 

de_ 
2tt' 



(BIO) 
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Now integrating over energy with the help of Eq. ( 4.12| ), we arrive at 



L a n (u;;V) = — ]T {2(p+l)g*(w)-M*(a;-2eV)-(p + l)(g + l)*(w + 2ey)}, 



p,g>0 

11 + 9 = 71-1 



L b n (ui;V) = — V (r + l)(s+l)[2#(u;)-#(cj-2eF)-#(w + 2eV)] 



(Bll) 



r,s>0 
r+3=n-2 



Evaluating the sums with the help of 



E 

p,q>0 

p+q=n— 1 

we get 



pq 



Y, (P + 1)5= E (r + !)(* + !) 



p,q>0 
p + q = n — 1 



t-,s>0 
r+s=n-2 



n + 1 



p,q>o 



Y: ( P +l)(g + l)=( n + 2 ), (B12) 



V) + L b n (u; V) = f^{4(« 2 - !)*(«) - (n - l)(2n - l)*(w - 2eV) - (n + l)(2n + + 2eV)}. (B13) 



Finally, substituting Eq. (BIS ) into Eq. (B7), we obtain 



S 2 



Analogously, 
6 2 



Tr K (e^Aoe-^Ao)" = i (2n 2 + 1) [tf - 2eV) + *(w + 2eV)l - 4(n 2 - l)*(w) ^ 

ip 2 =o 37r L J 

(B14) 

2 

-|-|(2n 2 + 1) - eV) + *(w + el/)] - 4(n 2 - l)*(w)}. 

(B15) 



¥=2=0 



APPENDIX C: SOLUTION OF THE USADEL 
EQUATION IN THE FIRST ORDER OVER A 



Below we solve the Usadel equation within the first order 
over Ad, neglecting, therefore, the effects of supercurrent. 



We will generalize here the method used in Ref. |l7| ] 
for the calculation of G a in the presence of interaction. 
Since we consider general case of an arbitrary interface 
transparency, the first step will be to find the interaction- 
induced correction to the "spectral angle" 9 E (r) that 
parametrizes [5]Jl7|l the retarded semiclassical Green func- 
tion in the N conductor, 



G R {E) = 



-le 



-up , 



ie 1 * sin( 

— COS0 



(CI) 



where the proximity angle 9 and the order parameter 
phase tp depend on the energy E and the space coordi- 
nate r. In N-S systems with a single superconductive ter- 
minal physical quantities do not depend on the phase on 
it, y(0), so below we put ip(0) = — 7r/2. In general, when 
interaction in the N conductor is present, normal current 
flowing through the N-S structure induces a supercurrent 
js oc V<y3, so that the phase ip(r) acquires some nontrivial 
distribution. However, this effect appears in the second 
order of expansion over Cooper interaction constant A^. 



1. Spectral angle 

The Usadel equation for the spectral angle 9e (r) has 
the form 



DW 2 9 E + 2iEsixi6 E + 2Acos0£; = 0, 



(C2) 



where the selfconsistency equation for the order param- 
eter reads 



A(r) 



Jo 



poo J£ 

-A,/ / dEtanh — Imsin6» E , (C3) 



and the interaction constant A<j at the energy scale tod is 
defined in Eq. (5.10). We assume that the main contri- 
bution to A comes from relatively high E, where E (r) 
is small everywhere. Hence, Eq. (02) can be linearized, 
and the last term can be neglected once we are interested 
in the first-order corrections over to Ga- The result 
is: 



5 (r) = A(E)K ((l-i)r/L E ), 



(C4) 
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where Le = WD/E, and A is determined by the bound- 
ary condition 



d6 E {r) 



dr 



Gt 



cos9e{t = d). 



(C5) 



The solution for A(E) is conveniently e xpres sed in terms 
of the function 9(f) defined in Eq. (fj~2l"l) : A(E) = 
Q(^B)/ha(-Lg/ d), w here t E — aln(uj d / E). Substituting 
9e{t) into Eq. (C3), using the identity / °° xK$(x)dx = 1 
and logarithmic slowness of A(tE) as a function of E, we 
find for A(r): 



A(r) 



A d L> 9(2aln(r/d)) 
r 2 In (r/d) 



(C6) 



The next step is to solve for 0E=o{r) with A(r) taken 
into account: 



d 2 9 1 dO _ 2\ d 9(2<zln(r/d)) 
dr 2 r dr r 2 In (r/d) 



cos 0. 



(C7) 



In terms of the variable £ = ln(r/d), this equation re- 
duces to 



d 2 9 9(2qQ 
_=2A d — ^— 



cos6> 



(C8) 



with the boundary conditions 0(£l) = and %(0) = 
—2a cos 0(0), where £ z = ln(L/d). To the lowest order in 
Xd one obtains 



9(0 = { 9(2< L ) - 2A d £ L 



1 + 2arysin9(2a£ L ) 9(2ary) 



l + 2a£ L sin9(2a£ L ) 
+ 2A d 



cos 



9(2a£i 



^ 9(2a?7) 
W~0 cos 



V_ 



9(2a£ z 



dry 



dry. 



(C9) 



Taking into account that 2o£l = a£ = t, and evalu- 
ating expression (C£) at £ = 0, we obtain the spectral 
angle near the interface, 9 d = 6(£ — 0): 



% = 9(t) 



A d C 



1 + <sin9(<) 
1 Q(xt) 



(1-x) cos [9(t)(l - a:)] dx. (CIO) 



2. Distribution function and Ga 

To find Ga we use the kinetic equation in the form 
developed in Ref. [[l7j. The current density near the N 
reservoir (at r = L) is: 



]n(R)=9 J dE 



dfi(E,r) 



dr 



(Cll) 



where f±(E,r) is the anomalous (bra nch imbalance) dis- 
tribution function, defined after Eq. ( |2.2| ) . Note that the 
Eq. ( |Cll| ) is valid only at the boundary with the normal 
lead where 9 = (otherwise it should be modified by 
the presence of the supercurrent j s oc A^). The function 
fi(E, r) obeys the equation 

VD(r, £7)V/i (E, r) = 2A(r) sin 9(E, r)f 1 (E, r). (C12) 

We will assume that the normal reservoir is biased by 
the voltage V, so that the whole distribution function is 
shifted: F(E,L) = tanh(£ + eVr z )/2T « tanh(£/2T) + 
t z (eV/2T) cosh" 2 (E/2T) and 



h (E,L) = eVc(E)= e ^ C o^r 2 ^. 



(C13) 



At T — > it is enough to consider Eq. (C12) at zero en- 
ergy, where D(0, r) = D, and the solution for 9(0, r) is 
given by 



9(0, r) = 9 d 



ln(L/r) 
ln(£/d)' 



(C14) 



with 9d given by Eq. (CIO). Then equation for fi(r) 
reduces to 



d£ 2 



-2A rf ^ sin 



(C15) 



We seek the solution of Eq. (C15) in the form 
/^(O + Ad/^O, where 



+ «d 1 



1 



(C16) 



The function c(-E) was defined in Eq. (CIS), and the 
parameter v d = V d /V should be determined from the 
boundary condition to the Usadel equation: 



that gives 



2a f i sin 9 d , 



1 



1 + 1 sin ( 



(C17) 



(C18) 



The general solution of Eq. (C15) for the function /j; 
can be written as 



(i) 



2G 



A (1) (0 = 



(C19) 



(i) 



9£ 



1 /" £i £i + 7?tsin6(t) 



L Jo 



l + *sin8(i) 



yfa)^. (C20) 



Calculating the total curre nt as 7jy = 2tt Lj^iL) where 



where y(£) = -2A d (6«)/0 sm[9(t)(l - ^l)]/i (0) (£)• 

Using Eq. (|C17| ) and the condition /^(fi) = 0, we find Jn{.L) is determined by Eq. ( |Cll| ) and introducing x — 
the constants Co and Ci, and obtain we obtain 

I 



j - 2n9 v 



t sin 8 r i 



4A d <| 



1 + tsin6» d 

Finally, the Andreev conductance is given by 



1 (1 + zt sin 6(f)) 2 0(xi) 
o (1 + < sin 6(f)) 2 arf - 



sin[6(i)(l - x)]dx 



i sin 9 r i 



G D l + tsvnO d 



1 (1 +xisine(t)) 2 Q(xt) 
o (1 + isin9(X)) 2 aT~ 



sin[0(t)(l - aj)]dx, 



(C21) 



(C22) 



with 9 d defined in Eq. ( |C1C| ). 

Eq. (C22) is an exact answer in the first order over 
Xd (in this approximation supercurrent does not y et mi x 
equations for /1 and /). In the t — > 00 limit of Eq. (C22), 
i. e., at Gt 3> Gd, the interaction- induced correction is 



Ga 
Gd 



- 2 i 1 L 
A d In — 

ir a 



(C23) 



and grows with the space scale, contrary to the results of 
Rcf. [O obtained for the ID geometry. 
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